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Abstr act. We prove the existence of knotted and linked thin vortex tubes for 
steady solutions to the incompressible Euler equation in R 3 . More precisely, 
given a finite collection of (possibly linked and knotted) disjoint thin tubes 
in R 3 , we show that they can be transformed with a C m -small diffcomorphism 
into a set of vortex tubes of a Beltrami field that tends to zero at infinity. 
The structure of the vortex lines in the tubes is extremely rich, presenting a 
positive-measure set of invariant tori and infinitely many periodic vortex lines. 
The problem of the existence of steady knotted vortex tubes can be traced 
back to Lord Kelvin. 
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1. Introduction 



The motion of particles in an ideal fluid in R 3 is described by its velocity field 
u(x,t), which satisfies the Euler equation 

d t u+ (u-V)u = -VP , divu = 

for some pressure function P(x,t). Equivalently, the field u satisfies the equation 

— dtu + u x uj = VB , div u = , 

where the field ui := curlu is the vorticity and B := P + ^\u\ 2 is the Bernoulli 
function. The trajectories (or integral curves) of the vorticity uj(x,t) for fixed t are 
usually called vortex lines. A solution u to the Euler equation is called steady when 
it does not depend on time. 

A domain in R 3 that is the union of vortex lines and whose boundary is an em- 
bedded torus is a (closed) vortex tube. The analysis of thin vortex tubes, in a sense 
to be specified below, for solutions to the Euler equation has attracted considerable 
attention. A long-standing problem in this direction is Lord Kelvin's conjecture |26] 
that knotted and linked thin vortex tubes can arise in steady solutions to the Eu- 
ler equation. This conjecture was motivated by results due to Helmholtz on the 
time-dependent case, which hinge on the mechanism of vorticity transport, and 
Maxwell's observations of what he called "water twists" . 

Kelvin's conjecture is basically a question on the existence of knotted invariant 
tori in steady solutions of the Euler equation. There is a considerable body of 
literature devoted to the analysis of topological and geometrical structures that 
appear in fluid flows, which has led to significant results e.g. on particle trajectories 
and vortex lines [2j [13j [20j [IT], on the relationship between the Euler equation 
and the group of volume-preserving diffeomorphisms [H [TUl HH [6], and on the 
connection of the helicity with the energy functional and the asymptotic linking 
number [3j HH [TT1 [27] . However, Kelvin's conjecture remains wide open, and indeed 
has been included as a major open problem in topological fluid mechanics in the 
surveys [251152]. 

There is strong numerical evidence of the existence of thin vortex tubes, both in 
the case of steady and time-dependent fluid flows. As a matter of fact, thin vortex 
tubes have long played a key role in the construction and numerical exploration of 
possible blow-up scenarios for the Euler equation, which in turn has led to rigorous 
results such as [9] . A particularly influential scenario in this direction is [24] , 
which discusses how an initial condition with a certain set of linked thin vortex 
tubes might lead to singularity formation in finite time. 

Our aim in this paper is to prove that there exist steady solutions to the Euler 
equation in R 3 having thin vortex tubes of any link and knot type. The steady 
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solutions we construct are Beltrami fields, that is, they satisfy the equation 

curl u = Xu 

in R 3 for some nonzero real constant A. As any Beltrami field satisfies the equation 
Au = — X 2 u in R 3 , it is apparent that the solutions we construct are real analytic 
but do not have finite energy (i.e., u is not in L 2 (R 3 )). However, our construction 
yields solutions with optimal decay at infinity in the class of Beltrami fields, which 
fall off as \u(x)\ < C/\x\. In particular, they are in L P (R 3 ) for all p > 3. 

The motivation to consider the class of Beltrami solutions to the Euler equation 
to address the existence of linked vortex tubes comes from Arnold's structure the- 
orem [21 Theorem II. 1.2]. Under mild technical assumptions, this theorem ensures 
that the vortex lines of a steady solution to the Euler equation whose velocity field 
is not everywhere collinear with its vorticity are nicely stacked in a rigid structure 
akin to those appearing in the study of integrable Hamiltonian systems. Heuristi- 
cally, this structure should somehow restrict the way the vortex lines are arranged; 
partial results in this direction have been shown in [12] . where it is proved that 
under appropriate (strong) hypotheses the vortex lines of steady solutions with 
non-collinear velocity and vorticity can only be of certain knot types. In contrast, 
Beltrami fields have been used recently to produce steady solutions of the Euler 
equation having a set of vortex lines diffeomorphic to any locally finite link |llj . 

Before stating the main result we need some definitions. We will say that a 
bounded domain of R 3 is a (closed) tube if its boundary is a smoothly embedded 
torus. Therefore, a vortex tube is a tube whose boundary is the union of vortex 
lines (or, equivalently, its boundary is an invariant torus of the vorticity field). A 
convenient way of constructing thin tubes is as metric neighborhoods of curves. 
Indeed, if 7 C R 3 is a closed curve, we will denote by 

(1.1) := {x e R 3 : dist(z,7) < e} . 

the tube of core 7 and thickness e. We are interested in the case where e is a small 
positive number, which corresponds to the case of thin tubes. Obviously any finite 
collection of disjoint tubes in R 3 can be isotoped to a collection of thin tubes of 
this form. 

Let us consider a finite collection of (possibly knotted and linked) disjoint thin 
tubes, constructed using metric neighborhoods of curves as above. Our main result 
in this paper is that, if the thickness of the tubes is small enough, this collection 
can be transformed by a C m -small diffeomorphism into a union of vortex tubes of 
a Beltrami field in R 3 . Furthermore, the structure of the vortex lines inside each 
vortex tube is extremely rich: firstly, the boundary of each vortex tube is far from 
being the only invariant torus of the flow, as the set of invariant tori has positive 
Lebesgue measure. Secondly, between any pair of these invariant tori there are 
infinitely many periodic vortex lines. Thirdly, there is a periodic vortex line which 
is close to the core of the initial tube and diffeomorphic to it. More precisely, we 
have the following statement, where (as always henceforth) all the diffeomorphisms 
are assumed smooth, and all curves and surfaces are smoothly embedded in R 3 . It 
should be emphasized that the thinness of the tubes is crucially used in the proof 
of the theorem. 

Theorem 1.1. Let 71,... ,Jn be N pairwise disjoint (possibly knotted and linked) 
closed curves in R 3 . For small enough e, one can transform the collection of pairwise 
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disjoint thin tubes 7^(71), . . . , 7^(7jv) by a, diffeomorphism $ ofM 3 , arbitrarily close 
to the identity in any C m norm, so that $[7^(71)], . . . , <&\J~<l{in)\ are vortex tubes of 
a Beltrami field u, which satisfies the equation curl it = Xu in R 3 for some nonzero 
constant A. 

Moreover, the field u decays at infinity as |Z? : 'u(a;)| < Cj/|a;| for all j and has 
the following properties in each vortex tube <£>[7^(7i)].' 

(i) In the interior o/$[7I(7i)] there are uncountably many nested tori invari- 
ant under the Beltrami field u. On each of these invariant tori, the field u 
is ergodic. 

(ii) The set of invariant tori has positive Lebesgue measure in a small neigh- 
borhood of the boundary 9$ [71 (74)]. 

(iii) In the region bounded by any pair of these invariant tori there are infinitely 
many closed vortex lines, not necessarily of the same knot type as the 
curve 74. 

(iv) The image of the core curve 73 under the diffeomorphism $ is a periodic 
vortex line of u. 

An important property of the structure of the vortex lines inside each vortex 
tube, as described above, is that it is stable in the following sense: on the one hand, 
it is robust under small perturbations of the field u, meaning that the trajectories 
of any field which is close enough to u in a sufficiently high C k norm have the 
same structure. On the other hand, the boundary of each vortex tube <9<I>[7^(7i)] 
is Lyapunov- stable under the flow of the Beltrami field u. It should be noticed too 
that from Property (iv) in Theorem 11.11 we recover the main theorem of [11] for 
finite links and improve it by ensuring that the solution falls off at infinity (while 
in [IT] we had no control at all on the growth of the solution at infinity). Besides, 
the proof of Theorem 11.11 shows that the vortex lines $(7^) are elliptic trajectories, 
and therefore linearly stable, while the vortex lines we constructed in |llj were 
hyperbolic, and thus unstable. 

After establishing his structure theorem, Arnold conjectured [2] that, contrary 
to what happens in the non-collinear case, Beltrami fields could present vortex lines 
of the same topological complexity as the trajectories of any divergence-free vector 
field. By KAM theory, typically these trajectories give rise to a set of invariant tori 
with positive measure and chaotic regions with homoclinic tangles between these 
tori. Theorem 1 1.1 1 is fully consistent with this picture, and proves the existence of 
the aforementioned positive-measure set of invariant tori. 

The paper is organized as follows. In Section [2] we will discuss the strategy of 
the proof of Theorem 11.11 This section, which serves as a guide to the paper, also 
allows us to explain the key difficulties that appear in the proof of this result but 
not in that of [llj . which require the introduction of new techniques and ideas, 
and make this paper considerably more involved. In Section [3] we introduce some 
objects associated with the geometry of a thin tube that will be used throughout 
the paper. In Section 3] we prove some estimates for the Laplacian in a thin tube 
with Neumann boundary conditions. These estimates are used in Section O to 
study harmonic fields in thin tubes. In Section [6] we construct Beltrami fields in 
thin tubes with prescribed harmonic part. These fields are analyzed further in 
Section!?] where we prove a KAM theorem for Beltrami fields in generic thin tubes. 
A Runge-type approximation theorem by global Beltrami fields tending to zero at 
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infinity is presented in Sectional With all these ingredients, in Section [9] we prove 
Theorem 11.11 The paper concludes with an easy application to the Navier-Stokcs 
equation, which we present in Section 1101 



2. Strategy of proof and guide to the paper 

To prove Theorem II .![ our basic goal is to establish the existence of a Beltrami 
field u, satisfying the equation 

curl u = Xu 

in R 3 for some nonzero constant A and falling off at infinity, such that the field u has 
a set of N invariant tori diffeomorphic to the surfaces {dT e (ji)}iLi} with 
being a set of prescribed (possibly knotted and linked) closed curves. We recall 
that Tei^i) is a metric neighborhood of the curve ji of small thickness e, as defined 
in Eq. By deforming them a little if necessary, we can assume without loss 

of generality that the curves ji are analytic. 

The basic idea behind the proof of Theorcm ll.ll is carried out in three interrelated 
stages. Firstly, we construct a local Beltrami field v, which satisfies the Beltrami 
equation in a neighborhood of each closed tube T e (ji) and has a set of invariant 
tori given by dl~ e (71), ■ • ■ , <97^(7_/v). Secondly, we prove that these invariant tori are 
"robust", meaning that suitably small perturbations of the local Beltrami field v 
still have a set of invariant tori diffeomorphic to dTe("fi), ■ ■ ■ , dTe (tjv)- To conclude, 
we show that the local Beltrami field v can be approximated in any C k norm by 
a global Beltrami field u, which satisfies the Beltrami equation in the whole space 
R 3 and falls off at infinity in an optimal way. The robustness of the invariant tori 
of the local Beltrami field ensure that u has a set of vortex tubes diffeomorphic to 
the initial configuration of thin tubes, and that in fact this diffeomorphism can be 
taken close to the identity in a C m norm. 

However, the implementation of this basic idea turns out to be extremely subtle. 
To understand why, one can start by noticing that the robustness of the invariant 
tori of the local Beltrami field relies on a KAM-type argument. To apply this 
KAM argument, we need a small perturbation parameter and some control on the 
dynamics of the local Beltrami field in a neighborhood of each invariant torus, 
which is required in order to show that the local Beltrami field is equivalent to a 
Diophantine rotation on the torus and satisfies a suitable nondegeneracy condition. 

To construct local Beltrami fields in a neighborhood of the tori dTd^i) with con- 
trolled behavior on these surfaces, it is natural to use some variant of the Cauchy- 
Kowalewski theorem for the curl operator with Cauchy data on the tori (e.g. [TTJ 
Theorem 3.1]). This would lead to a local Beltrami field defined in a small neigh- 
borhood Hi of each torus dT t {^i). Unfortunately, the Runge-type approximation 
theorem we prove in this paper does not allow us to approximate a local Beltrami 
field defined in Vt\ U ■ • • U f2jv by a global Beltrami field, since the complement 
R 3 \(r2i U ■ ■ ■ U Qn) has compact connected components. As is well known, this 
is not just a technical issue, but a fundamental obstruction in any Runge-type 
theorem. 

Therefore, to construct the local Beltrami field v we will not consider a Cauchy 
problem but a boundary value problem for the curl operator in each tube. This 
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leads to a vector field that satisfies the Beltrami equation 

curlv = Aw 

in a neighborhood of the tubes (because the boundaries are analytic), which allows 
us to apply our Runge-type approximation theorem. In this boundary value prob- 
lem one can prescribe the normal component of v at each boundary dTeirii), so that 
by setting it to zero we can ensure that each boundary is an invariant torus of the 
local Beltrami field. However, we have no control on the tangential component of 
v on each torus, so in principle the local Beltrami field does not necessarily satisfy 
the dynamical conditions our KAM-type theorem requires for the invariant tori to 
be preserved under small perturbations. 

To overcome this difficulty we resort to a careful analysis of harmonic fields in 
thin tubes. Indeed, we show that in the above boundary value problem for the 
local Beltrami field, one can also prescribe its harmonic part (that is, the L 2 pro- 
jection of v into the space of harmonic fields with tangency boundary conditions), 
and for small A one expects the local Beltrami field to behave essentially as if it 
were harmonic. Hence, we exploit the small parameter e (that is, the thickness of 
the tubes) to extract detailed analytic information about harmonic fields through 
energy estimates, and then utilize this knowledge to compute the dynamical prop- 
erties of the local Beltrami field that are required in the KAM-type theorem under 
the assumption that the Beltrami parameter A is suitably small (in fact, of order 
e 3 ). A key ingredient in the computation of these dynamical properties is a set of 
energy estimates for the local Beltrami fields and all its derivatives that are optimal 
with respect to the geometry of the tubes, i.e., the parameter e. 

The local Beltrami field can now be approximated in any C k norm by a global 
Beltrami field u that falls off at infinity as C/\x\ using a Runge-type theorem. When 
using the KAM argument to guarantee that u still has a set of invariant tori diffco- 
morphic to {dTeiji)}, one has to face the problem that the local Beltrami field is in 
fact a small perturbation of a field that does not satisfy the nondegeneracy condi- 
tion (equivalent to a rotation of the disk with constant frequency), which requires 
a fine analysis of the terms controlled by each small parameter: the thickness e, the 
Beltrami constant A and the error in the Runge-type approximation. 

We shall next present a sketch of the proof of Theorem 11.11 where we explain 
the different intermediate results that are needed in the demonstration and their 
interrelations. This short sketch is also intended to serve as a guide to the paper. 
For the sake of clarity, we will divide the proof in three stages: 

Stage 1: Construction of the local Beltrami field. The local Beltrami field 
v is obtained as the unique solution to certain boundary value problem for the 
Beltrami equation. Our goal is to estimate various analytic properties of this field, 
and for this it is natural to introduce coordinates adapted to a Frenet frame in each 
tube 71 = Teipfi), which essentially correspond to an arc- length parametrization 
of the curve 74 and to rectangular coordinates in a transverse section of the tube. 
Thus we consider an angular coordinate a, taking values in E>\ := R/£Z (with I the 
length of the curve 7,), and rectangular coordinates y = (2/1,2/2) taking values in 
the unit 2-disk B 2 . Details are given in Section^ 
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To extract information about the local Beltrami field, a useful tool is the bound- 
ary value problem for the Laplacian on scalar functions with zero mean and zero 
Neumann boundary conditions in the thin tube T e : 

Ai/j = p in T e , d v ip\ar t = , J ip dx = . 

When written in the natural coordinates (a, y), we obtain a boundary value prob- 
lem in the domain §^ x D 2 , the coefficients of the Laplacian in these coordinates 
depending on the geometry of the tube strongly through its thickness e and the 
curvature and torsion of the core curve ji . 

It is clear that we have the H k estimate 

IMIff*(T e ) < C e,k\\p\\Hk-2(T e ) ■ 
However this estimate, where the constant C e .fe depends on e in an undetermined 
way, is far from being enough to compute the dynamical quantities for the local 
Beltrami field that are needed in the KAM-type theorem. Therefore, in Section^ 
we prove energy estimates for the function tp and its derivatives that are optimal 
with respect to the parameter e. In particular, in order to be able to compute 
the desired dynamical quantities for the local Beltrami field later on, it is crucial 
to distinguish between estimates for derivatives of ip with respect to the "slow" 
variable a and the "fast" variable y. The estimates for the function ip and its 
derivatives that we will need are stated in Theorem 14. f II 

These estimates are immediately put to work in Section[5\to derive a perturbative 
expression for the harmonic field in each thin tube 71 and small e. (Of course, 
the problem degenerates for e = and the asymptotic results we prove do not 
correspond to a Taylor expansion.) We will use the notation h for the harmonic 
field in 71, which is unique up to a multiplicative constant. In fact, we need to 
compute h up to corrections that are suitably small for small e; as before, different 
powers of e are required for the slow and fast components of the field. To first 
order, the harmonic field can be written in the coordinates (a, y) as 

(2.1) h = d a +r(a)(y 1 d 2 ~y 2 d 1 )+0(e), 

where r stands for the torsion of each curve ji and 0(e) denotes some vector 
field whose components in these coordinates are bounded by a multiple of e in the 
C k (S>] x D 2 ) norm. By Eq. (|2.1I) , the field h is an 0(e) perturbation of the rotation 

of constant frequency given by the total torsion J Q r(a) da. From the point of view 
of KAM theory, this is a very degenerate case, so it is not hard to guess that one 
needs to compute h (at least) up to order 0(e 3 ), as we do in Theorem 15.11 As a 
matter of fact, we will see later on that the nondegeneracy condition that appears 
in the KAM theorem is that a certain quantity, called the normal torsion, must be 
nonzero, and that it actually vanishes modulo terms of order 0(e 2 ), which justifies 
why we need good estimates for h. 

In Section&we show that, for any nonzero constant A that is small enough (say, 
smaller than some fixed, e-independent constant), there is a unique vector field that 
is tangent to the boundary of each tube T e , satisfies the equation 

(2.2) curlv = Xv 

in the tube and whose harmonic part (i.e., the L 2 projection on the space of har- 
monic fields) is h. (This is, of course, different to what happens in the case of 
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compact manifolds without boundary, where Beltrami fields are always orthogonal 
to harmonic fields.) We also prove estimates that measure how the field v becomes 
close to h in the C k norm for small e and A. The key result is Theorem 16.81 which 
will be crucial in verifying the conditions in the KAM argument for the preservation 
of invariant tori. 

Stage 2: Preservation of the invariant tori. In Stage 1, for small enough e 
and A, we have constructed a local Beltrami field v, which satisfies Eq. (|2.2|) in 
a neighborhood of \J i=1 T e (ji) and has a set of invariant tori given by dT e (^i)- 
Furthermore, the estimates we have proved provide a very convenient expression 
for the local Beltrami field, up to terms that are of order e 3 (when A = e 3 ) in a 
suitable sense. 

In Section [?] we analyze the robustness of the invariant tori using the Poincare 
map n defined by the local Beltrami field at a transverse section of each tube. In 
each tube, II is then a diffeomorphism of the disk that preserves certain measure, 
while the intersection of each invariant torus with the transverse section is an in- 
variant circle of the Poincare map. The persistence of the invariant tori will rely 
on a KAM theorem for the Poincare map (Theorem IT.6|) that applies to individual 
invariant circles. 

To apply the KAM theorem, one has to verify that the rotation number of II on 
the invariant circle is Diophantine and that II satisfies a nondegeneracy condition 
(namely, that its normal torsion is nonzero) . Computing these dynamical quantities 
using the estimates for the local Beltrami field is nontrivial. On the one hand, the 
rotation number depends on the behavior of the trajectories for arbitrarily large 
times, so a delicate treatment is required in order to control uniformly the effect of 
the C(e 2 ) terms in the vector field. On the other hand, the normal torsion, which 
is required to be nonzero, turns out to be zero modulo C(e 2 ), since the Poincare 
map is a small perturbation of a constant frequency rotation of the disk, which is 
a highly degenerate case from the point of view of KAM theory. 

The estimates proved at Stage 1 are tailored to permit us to overcome these diffi- 
culties, yielding expressions for the rotation number wri and the normal torsion Afu 
that depend on the geometry of each curve 7$ through its curvature k and torsion r 
(see Theorems El and [7?HJ): 



These expressions allow us to prove the main result in this Stage 2, which is that 
for "generic" curves 7, the hypotheses of the aforementioned KAM theorem are 
satisfied, so that the invariant tori of the local Beltrami field v are robust: if u is a 
divergence-free vector field in a neighborhood of the tubes that is close enough to v 
in a suitable sense (e.g., in a high enough C k norm), then u also has an invariant 
tube diffcomorphic to each Te("fi), and moreover the corresponding diffcomorphisms 
can be taken close to the identity. This is proved in Theorem 17. 101 
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Stage 3: Approximation by a global Beltrami field. In Section\^we prove a 
Runge-type approximation theorem for Beltrami fields that decay at infinity (The- 
orem |531). More precisely, we will show that the local Beltrami field v, considered 
in the previous stages and defined in a neighborhood of the thin tubes 71 (74), can 
be approximated in any C k norm by a Beltrami field u that falls off at infinity as 
\D^u{x)\ < Cj/\x\. This decay is optimal for the class of Beltrami fields. 

The proof of this theorem consists of two steps. In the first step we use functional- 
analytic methods to approximate the field field v by an auxiliary vector field w that 
satisfies the elliptic equation Aw = —\ 2 w in a large ball of M 3 that contains all 
the tubes. In the second step, we define the approximating global Beltrami field u 
in terms of a truncation of a suitable series representation of the field w, ensuring 
that u has the desired fall-off at infinity. 

Completion of the proof of the main theorem. In Section [5| we use all 

previous results to complete the proof of Theorem ll.il Indeed, from the robustness 
of the invariant tori of the local Beltrami field, it immediately follows that the 
global Beltrami field u has a set of thin vortex tubes equivalent through a C m - 
small diffcomorphism $ to {7^(7i)}, provided e is small enough. More precisely, 
what is proved is that there is some constant eo, which depends only (albeit in 
a rather nontrivial way) on the geometry of the curves 7$ and on the allowed 
smallness for the diffcomorphism $, so that the statement of the theorem holds for 
any thickness e ^ eo. 

The remaining properties of the vortex lines stated in Theorem ll.ll are established 
in this section too, using results derived throughout the paper. In particular, we 
show that near the curve % there is an elliptic periodic trajectory of the field u. 
As a side remark, notice that this elliptic periodic trajectory is obviously linearly 
stable, but is not granted a priori to be Lyapunov stable. Therefore, as shown by 
the counterexample of Anosov and Katok (cf. e.g. [E]), it is not guaranteed that 
there are invariant tori of the field in a neighborhood of the elliptic trajectory. A 
careful (and nontrivial) analysis of the dynamics near the elliptic trajectory for 
small e and A would be required to prove the existence of these tori. 

Before beginning with the technical part of the paper, it is convenient to provide 
a short comparison between the proof of Theorem 1 1 . 1 1 and that of the main theorem 
in Ref . [IT] , where we showed that there are steady solutions to the Euler equation 
with a set of vortex lines diffeomorphic to any given link. In this reference, the 
proof was also based on the construction of a local Beltrami field with a "robust" 
set of periodic trajectories, which was then approximated by a global Beltrami field. 
However, the implementation of this basic principle is totally different. 

To begin with, the robustness of periodic trajectories in [UJ relies on the hy- 
perbolic permanence theorem (which is essentially an application of the implicit 
function theorem) instead of a considerably more sophisticated KAM argument. 
To construct a local Beltrami field with prescribed periodic trajectories 7* that are 
hyperbolic, as required by the permanence theorem, it is enough to prove a suitable 
analog of the Cauchy-Kowalewski theorem for the curl operator: indeed, since the 
field is divergence-free, it is enough to assume that the field used as Cauchy datum 
is exponentially contracting into the curve 7$ on the Cauchy surface to obtain a 
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local Beltrami field that has ji as a hyperbolic periodic trajectory. All this is in 
strong contrast with the case of vortex tubes, where the construction of the lo- 
cal solution with a robust set of prescribed invariant tori requires the analysis of 
the boundary value problem and the KAM argument described in Stages 1 and 2 
(Sections H to [3). 

The approximation theorem we use in this paper is also different from the one 
employed in [IT] . The reason for this is that in [IT] we had no control on the growth 
of the global Beltrami field at infinity, even in the case of connected links. On the 
contrary, the approximation theorem we prove in this paper (Theorem 18. 3p yields 
Beltrami fields with optimal fall-off at infinity. The proof of these approximation 
theorems is considerably different: the old theorem is based on an iterative scheme 
that uses a theorem by Lax and Malgrange and works for the curl operator in any 
Riemannian 3-manifold, while the new one is based on different principles and takes 
advantage of the geometry of Euclidean 3-space to ensure that the global Beltrami 
field falls off at infinity. 

Notice that the main theorem in [11] applies to locally finite sets of curves, while 
in Theorem 1 1.1 1 in this paper we can only take finite sets of curves 71, ... , 7 at. Some 
comments in this regard are made in Remark 19.21 

3. Geometry of thin tubes 

In this section we introduce some notation, including a coordinate system, that 
will be used throughout the paper to describe functions and vector fields defined 
on thin tubes. These tubes will be characterized in terms of the curve that sits on 
its core and its thickness e, which is a parameter that will be everywhere assumed 
to be suitably small. 

Let us start with a closed analytic curve with an arc-length parametrization 
7 : Sj — > M 3 , with Sj := R/£Z (throughout the paper, when the period is 27T we 
will simply write S 1 = S^tt)- This amounts to saying that the tangent field 7 has 
unit norm and I is the length of the curve. We will abuse the notation and denote 
also by 7 the curve parametrized by the above map (i.e., the image set 7(8^)) . 

Let us denote by T € = a metric neighborhood with thickness e of the 

curve 7, that is, 

T e := {iel 3 : dist(z,7) < e} . 

This is a thin tube having the curve 7 as its core. It is standard that, for small e, the 
boundary dl~ e is analytic. The normal bundle of the curve 7 being trivial |21) . one 
can associate to each a£§J two orthogonal unit vectors ej(a) in M 3 perpendicular 
to the curve at the point 7(a). For convenience, we will make the assumption that 
the curvature of the curve 7 does not vanish, which allows us to take ei(ct), 62(a) 
as the normal and binormal vectors at the point 7(a). It is well known that the 
assumption that the curvature does not vanish is satisfied for generic curves [5] 
p. 184] (roughly speaking, "generic" refers to an open and dense set, with respect 
to a reasonable C k topology, in the space of smooth curves in R 3 ). 

Using the vector fields &j[oi) and denoting by ID 2 the two-dimensional unit disk, 
we can introduce analytic coordinates (a,y) G x B 2 in the tube 71 via the 
diffeoniorphism 

(a, y) 1 v 7(a) + eyiei(a) + ey 2 e 2 (a) . 
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In the coordinates (a, y), a short computation using the Frenet formulas shows that 
the Euclidean metric in the tube reads as 



(3.1) ds 2 = Ada 2 + 2e 2 T(y 2 dyi - yi dy 2 ) da + e 2 (dy 2 + dy 2 ) , 

where k = k(o) and r = r(a) respectively denote the curvature and torsion of the 
curve, 

(3.2) A := (1 - en Vl ) 2 + (er) 2 \y\ 2 , 



and \y\ stands for the Euclidean norm of y = (2/1,2/2)- As is customary, we will 
denote by g%j and g lJ the components of the metric tensor and its inverse, respec- 
tively. 

We will sometimes take polar coordinates r £ (0, 1), 9 £ S 1 := R/27rZ in the 
disk B 2 , which are defined so that 

yi = r cos 6 , y 2 = r sin 9 . 

The metric then reads 

(3.3) ds 2 = Ada 2 ~ 2e 2 rr 2 d9 da + e 2 dr 2 + e 2 r 2 d9 2 , 

where we, with a slight abuse of notation, still call A the expression of (|3.2| in 
these coordinates, i.e., 

(3.4) A := (1 - enr cos 9) 2 + (err) 2 . 

Notice that the coordinate r is simply the distance to the curve 7. For future 
reference, we record that the volume measure is written in these coordinates as 

(3.5) dV := Bdady = Br da dr d6 
up to a factor of e 2 , where dy := dyi dy 2 and 

(3.6) B := 1 — enyi = 1 — enr cos 9 . 

4. Estimates for the Neumann Laplacian in thin tubes 

In this section we will derive some estimates for the Laplace equation Ai/> = p in 
the thin tube T e with zero Neumann boundary conditions. As we will see, to gain 
control on the function ip it is convenient to attack this equation in the coordinates 
(a, y) £ E>j x D 2 that we introduced in Section [3j in terms of which the Laplace 
equation can be written as 

(4.1a) AV> = p in §JxD 2 , d v ip = on S,} x <9D 2 , 

where the Laplacian A is now interpreted as an e-dependent differential operator 
in the variables (a,y) or {a,r, 9). In order to ensure the existence of solutions to 
this equation, we suppose that J pdV = 0, which allows us to uniquely determine 
the solution ip by demanding that it also has zero mean: 

(4.1b) J i;dady = 0. 

(Here and in what follows, we omit the domain of integration when it is the whole 
domain S} x D 2 . We could have used the measure dV in the above formula too, 
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but this choice is slightly more convenient.) For future reference, we record here 
the expression of A in the variables (a, r, 8): 

1/ 1, A , \ 1, 2t, t' - erUr' - k't) cos 8 , 
= [i>„ + -A + ^B2^ee ) + -^Vw + -p Vafl + ^3 

(4-2) 

1 (K sin 6(B 2 — (err) 2 ) kcosO \ er{n' cos 8 — tk smd) 

+ 7 \ 7& ^ e - -b - ^ J + P ^ • 

As usual, we denote partial derivatives by subscripts when there is not risk of 
confusion. 

Given a subset ft C §>] x B 2 , we will use the notation 



|n := ( / "0 2 dady 
n 



1/2 



for the L 2 (f2)-norm of the function ip, omitting the subscript when Q is the whole 
domain S] x D 2 . In this section we will use the notation 

(4-3) H\\l } := \\d^\\ 2 + ^t 

for a homogeneous Sobolev norm in which the derivatives associated with the "small 
directions" of the thin tube are weighted with an appropriate e-dependent factor. 
The usual H k norm of the function ijj(a,y) will be denoted by 1 1 1 1 - 

In what follows, we will assume that ip is a solution of the Laplace equation (|4.ip 
with zero mean and zero Neumann boundary conditions, which ensures that 

(4.4) J g ij diip djip dV = - J pipdV 

for any ip £ ff 1 (§] x D 2 ). As usual, the scripts i, j range over the set of coordinates 
{a, 7/1, j/2}, and summation over repeated indices is understood. It is important to 
notice that, for any function ip, 

(4.5) J g i Jd iV >d j <pdV=(l + 0(e)) |M|^ , 

so that the Hi norm is essentially a more convenient way of dealing with the natural 
H 1 norm associated with the metric. We will use this identity many times in this 
section without further comment. 

The structure of this section is the following. We will start by estimating the 
I? norm of the derivatives of the function ip with respect to the "slow" variable a 
(Subsection ^. ip . The proof of these estimates is standard. We recall that the reason 
why a is called the slow variable is that it parametrizes the "large" direction of the 
thin tube T t , as opposed to the "fast" variable y, which is obtained by rescaling 
the small section of the tube. Estimates for the derivatives of i/j with respect to the 
"fast" variable y with optimal dependence in the small parameter e are presented 
in Subsection 14.21 They are more complicated to obtain, basically because one has 
to consider an auxiliary function in order to get rid of the contributions to tjj that 
only depend on the slow variable. The resulting estimates for ip, which we will 
often use in forthcoming sections, will be stated in Subsection 14.31 
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4.1. Estimates for derivatives with respect to the "slow" variable. In this 
subsection we will prove H k estimates for the derivatives of the function t/j with 
respect to the "slow" variable a (cf. Proposition I4.3[) . We will begin with the 
following proposition, where we estimate the L? and H} norms of the function 
ip. As is customary, throughout this article we will use the letter C to denote 
e-independent constants that may vary from line to line. 

Proposition 4.1. The function ip satisfies the estimate 

\m+u\\Ai<c\\p\\. 

Proof. By the expression of the metric in the coordinates (a, y), it is clear that for 
small enough e one has 

g^d^d^dV^ {l + 0(e)) J (v£ + ^L)dad W >i J ^l + \D y ^\ 2 ) dady . 

The rightmost term in the inequality is bounded from below by ^/ii||'0|| 2 , where pi 
stands for the first nonzero Neumann eigenvalue of the flat solid torus: 

Mi := inf | J (<f£+\ D v<P\*) dad V : V e C°°(SjxB 2 ), J ^dady = 0, \\(p\\ = 1 j > . 
Hence we infer that 

(4-6) M < C'U\\^ 

for some C > 0. 

To conclude, we can now use the weak formulation of the equation (|4.4I) with 



<p = ip and the Cauchy-Schwartz inequality to derive that 

< (l + 0(e)) J g ij d^d^dV^ (l + O^WpWUW^CWpWM^. 

The proposition then follows from this inequality and the estimate (|4.6[) . □ 

In the following lemma we record an elementary inequality that will be of use 
several times in this section. 

Lemma 4.2. Let D a y denote the tensor of I th order derivatives with respect to the 
variables (a, y). For any well-behaved (e.g., smooth) functions ip, <p, x, x on S\ xD 2 , 
one has 



D l a jBg^)d^d^xx\ dady sC eQ U g*> d i9 d j9X dV\ (J g^d^djPxdV 
In particular, 

J \D l aty (Bg ij )d i (pd j (p\dady Cie\\^\\jjx H^Hhi ■ 

Proof. It is an immediate consequence of the expression for the metric and the 
function B in the coordinates (a,y) (see Eq. (|3.1|1 ^ and the Cauchy-Schwartz in- 
equality. □ 

We are now ready to prove the estimates for the derivatives of ip with respect to 
the slow variable a that we will need in this paper: 
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Proposition 4.3. The k th partial derivative of the function ip with respect to the 
angle a satisfies 

ii^+vh^ <c k \\ P \\ H >, 

where k is any nonnegative integer. 

Proof. Let us take ip — ip aa in Eq. (|4.4j) and integrate by parts to get 

g ij diip a djip a dV = J pip aa dV - J d a {Bg 13 ) diip djip a da dy . 

The LHS is bounded from below by (1 + 0(e)) H^ccH^L , while from the definition of 
the norm H} and Lemma l4~2l it stems that 



pip aa dV 



<(l + 0(e))||p||||^ a ||^, 



d a {Bg lJ ) diip djip a da dy 
Using Proposition 14. II to estimate HV'lliji in terms of \\p\\, we infer that 

ll^l&i ^c\\p\\\\i> a \\ Ai , 

which readily implies the desired bound for k = 0. When A; is a positive integer, the 
proof is totally analogous and can be obtained by induction on k using (p = da k+2 ip, 
the only difference being that one needs to estimate the term J pip as 

k+l 

pdl k+2 ^dV ^CWpWh^W^Whi ^C-IIpII^+CIIpII^H^IU; 
J=l 

by the induction hypothesis. 



□ 



4.2. Estimates for the "fast" variables. To estimate the derivates with respect 
to the "fast" variable y in an optimal way, it is crucial to ensure that the terms that 
only depend on a are not considered when estimating the norms of the function. 
A convenient way of doing this is by considering the auxiliary function 

(4.7) i>{oL,y) :=ip(a,y) / tp(a,y')dy' , 



which is obtained from ip by subtracting its average in the fast variable. It should 
be emphasized that, if we were not to subtract this average, the estimates we would 
obtain would not be strong enough for our needs in later sections. The key estimate 
in this subsection is Theorem [ 



An immediate observation is that, of course 
(4.8) Di y d& = D{d k J 

whenever the number j of derivatives we take with respect to y is greater than zero, 
and that ip has zero mean: 

J ipdady = . 

Moreover, the function ip satisfies the equation 

AiP = p in S} x © 2 , d u ip = on S} x <91D> 2 , 
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with 



p{oi, y) ■= p(a, y)--A [ ip(a, y') dy' ■ 



The estimates for ip we derive in the previous section guarantee that the norm 
of p can be bounded by a multiple of the norm of the initial source term p: 

Proposition 4.4. The H norm of p is bounded by 

WPWh* < Ck\\p\\m ■ 

Proof. Observe that the action of the Laplacian (which we compute in the coordi- 
nates (a, j/)) on the function ip(a,y') is 

A</<a, y') = (i + 0(e)) afcKct, y') + 0{e) 9 Q v(«, y') , 

where C(e n ) here stands for an e-dependent quantity Q(ot, y) bounded as 

HOIIh* ^C k e n 
for all k. Therefore, one finds that 



\D j y d k a p\\ = 



Dld k aP --Dld k a A [ iP(a,y')dy' 
n Jo 2 

D 3 y d k a p-- [ d k a + 2 i>(a,y')dy'+J2 



0(e)d l a ip(a,y>)dy' 



< \\D y d k p\\ + (± + ffe) \\d k+2 yj\\ + Ce J2 O 
<C\\p\\ Hk , 

where in the last step we have used Propositions 14.11 and 
follows. 

To derive energy estimates for ip, we will use that one obviously has 



The claim then 
□ 



(4.9) 



pip dV 



for all ip £ x D 2 ). Our first result will be an estimate for the L 2 norm 

of ip and Dyip. While we can readily derive a bound for these quantities using 
Proposition 14. 1[ the estimates we prove here are much sharper for small e. This 
will be crucial in the derivation of optimal estimates for ip. 

Proposition 4.5. The function ip satisfies the H 1 estimates 



m + \\D y ^\\^Ce 2 \\p\\ 



\d a $\\ < Ce|MI 



Proof. Choosing ip — ip in Eq. (|4.9[) and in view of the expression of the coefficients 
of the metric (|3.1[) . one immediately obtains that 



(4.10) 



H} 



;i + 0(e))||p||||^|| 



Since tp(cx, •) has zero mean in the disk B 2 for any fixed a, by Poincare's inequality 
there is a positive constant C, independent of e and a (namely, the first nonzero 
Neumann eigenvalue of the disk), such that 

\D y tp(a,y)\ 2 dy > C ip(a,y) 2 dy. 
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Integrating this inequality in a, the H 1 norm of ip can be estimated as 
Wi>Wh>^2 i \D y 4>{a,y)\ 2 dyda^- 2 U\\ 2 - 



Together with Eq. (l4~TU|) . this yields \\ip\\ Ce 2 \\p\\ and H^H^i Ce\\p\\, so the 
claim follows directly from Proposition 14.41 □ 



It is particularly easy to derive preliminary estimates (which will be instrumental 
in the proof of Proposition HIT]) for the derivatives of tp with respect to a thanks to 
Proposition ^. 31 Notice these bounds will be substantially improved later on. 

Proposition 4.6. The derivative of ip with respect to a satisfies 

\\D y d k a +1 n < c 4p\\h* , \\d k a +2 4>\\ < c k \\ P \\ Hk 

for any nonnegative integer k. 

Proof. The claim is an immediate consequence of the definition of xp (cf. Eq. (|4.7|l ) 
and Proposition I4.3[ since obviously the L 2 (S] x ED 2 ) norm of the function 

di [ iP(a,y')dy' 



is bounded from above by C||9^||, which was in turn estimated in the aforemen- 
tioned proposition. □ 

We are ready to show that the second derivatives of ip with respect to the fast 
variable y are bounded by a factor of order e 2 . The proof of these estimates makes 
essential use of Propositions 14.51 and 14.61 

Proposition 4.7. The second derivatives of the function ip with respect to the fast 
coordinates are bounded by 

\\D 2 y M\ < Ce 2 \\p\\ . 

Proof. We will denote by the two-dimensional disk of radius R, R being a fixed 
real smaller than 1. Let us start by proving interior estimates. For this, we will 
denote by d a the derivative along a y-direction (that is, d\ or <9 2 ), and consider a 
smooth function x(M) equal to 1 for \y\ < R and equal to zero in a neighborhood 
of <9B 2 . Taking ip = d a (x 2 d a ip) in Eq. (|4.9[) (here and in what follows we will not 
sum over the index a) and integrating by parts, one readily obtains 

(4.11) I 2 := J g ij did a i> djdai) x 2 dV = h-I 2 -h, 

where 

h := [ pd a (x 2 d a 4>)dV , h ■= I d a (BgV)dSd 3 (x 2 d a 4>)dady 

2\ 



h ■= J da^g^djix^d^dV . 

Let us estimate these integrals. The first one can be easily controlled using the 
Cauchy-Schwartz inequality and Propositions 14.41 and 14.51 



Px'd^dV 



pd a ( X 2 )d a ^dV 



^C\\p\\(\\ X d 2 J\\ + \\d a ^\\ 



< Ce\\p\\I + Ce 2 \\p\\ 
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The integral I2 can be controlled using an analogous argument, Lemma 14.21 and 
Jensen's inequality. This leads to the estimate 



I/2I 



d a {Bg lJ )x diip djd a ip da dy 

^ Ce\$\\ Al I + C\\d a j>\\ C f (d a (Bg^)dSd jX y 'dady 

da dy 



d a {Bg ij )dii, dj^daj dady 



^Ce'\\p\\I+Ce 2 \\p\\ / d a (Bg^)d^d jX 
< Ce 2 \\p\\I + C<?\\p\\\$\\ 6i \\xWhi < Ce 2 \\p\\I + Ce 3 \\p\\ 2 , 

where in the fourth line we have used that Hxlliji — li-Dj/Xll/ 6 = C/e. A similar 
argument shows that 

\h\^C\\da$\\\\xUiI<Ce\\p\\I. 
Feeding these bounds into Eq. (|4.11|) . we obtain 

/ 2 ^C £ ||p||/ + Ce 2 ||p|| 2 , 

so that I ^ Ce||p||. From the definition of the function x an d the identity (|4. 5[) it 
then follows that 

(4-12) ll^llsjxDj <Ce 2 ||p||, 

which is the desired interior estimate. 

To prove the estimates up to the boundary, we begin by showing that the 
Hi norm of 

d e ip = yi d 2 ^ - V2 dxip 
is bounded in terms of ||p||. For this, we will find it convenient to take a smooth 
function x(M)i equal to 1 for \y\ > R and vanishing in a neighborhood of the 
origin, and use polar coordinates throughout without further notice. If we now 
take <p = X 2 4 > 60 in Eq. (|4.4[) and integrate by parts, we readily find that 



g^ dSe drfe X 2 dV = / pi>ee X 2 dV ~ d e (Bg* ) d$ d 3 ( X 2 ^) da dy 



where the indices i,j now range over the set {r, a}. Notice that the reason we 
are now using polar coordinates is that dg does not commute with the derivatives 
with respect to (2/1,2/2)- Arguing as above, one finds that 

t g lJ dS e d^ e x 2 dV ^ Ce\\p\\ 

which ensures that 

(4-13) H&flllsJxA* + Hre\y tXAR < Ce 2 \\p\\ , 

where A# := J3 2 \p 2 R is the annulus of inner radius R. 

To estimate the derivative i£> rr , it now suffices to isolate this quantity in the 
equation At/j = p. From the expression of the Laplacian in these coordinates (|4.2[) 
it is apparent that for r > R one can write ip rr as 

e - 2 Vv r -p = 0{e- 2 ) (\$ee\ + + Oie- 1 ) $ e + 0(1) (\4> ae \ + $ aa \ + \$ a \) ■ 
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From Propositions I4.5H4.7I and the estimates (|4. 1 3[) , it then follows that 

IIVVrllsJxAH < Ce 2 \\p\\ . 

This yields the desired boundary estimates, thus completing the proof of the propo- 
sition. □ 

The results we have established so far show that the derivatives of -tp can be 
bounded in terms of the source p as 

(4.14a) |H| + ||^|| + ||D^||<Ce 2 ||p||, 

(4.14b) \\d a i>\\ + \\D v d a $\\ sc Ce\\p\\ , 

(4.14c) \\d 2 J\\ sc C\\p\\ . 

However, having in mind applications in forthcoming sections, we would rather have 
estimates where the RHS always has a factor of e 2 . 

In the following theorem we show that this can be achieved by replacing the 
L 2 norm of p by its H 1 norm (thus using estimates that are weaker in terms of the 
gain of derivatives) , and provide a generalization for higher derivatives. It is worth 
emphasizing that both the estimates (|4.14l) and those in the following theorem are 
optimal with respect to its dependence on the small parameter e, as can be checked 
easily. 

Theorem 4.8. For any nonnegative integers j and k we have the bound 

\\Did*M^C 3k e 2 \\p\\ HJ+kl 

where J := max{j — 2, 0}. 

Proof. The proof proceeds by induction on J + k, that is, on the number of deriva- 
tives in the RHS of the inequality. The case J + k = follows from Propositions 14. 51 
and [473 To avoid cumbersome notations that might obscure the argument, we will 
sketch the procedure for the case J + k = 1, where one has to estimate d a ^, D y d a ^ 
(improving the bounds in Propositions 14.51 and I4.6[) , Dyda^p and Dy'tp. Once this 
case has been worked out in detail, it is straightforward to prove the general result 
using an induction argument. 

Let us begin by estimating the quantities having derivatives with respect to a, 
that is, 9 Q i/>, Dyd a 'ip and D 2 d a ip. An easy computation using the expression of the 
Laplacian (|4. 2|) shows that the commutator 

q := A(d a i>) - d a (Aij) 

can be written as 

q = 0{e- 1 ) Dl^P + 0(e) + 0(1) D v d a ^ + O^Dy^ + 0(e)j a . 

Now from the H 2 estimates proved in Propositions I4.5H4.71 one obtains that the 
L 2 norm of the commutator is bounded by 

y|<Ce|H|. 

The function d a tp obviously satisfies the zero-mean condition 

(4.15a) / d a ip(a,y')dy' = 0, 

jd 2 
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the boundary condition 

(4.15b) d u (d a ip~)=0 

on B>\ x <9ID 2 and the equation 

(4.15c) A(d a ip) = d a p+g. 

As the L 2 norm of the RHS is bounded by C||/9||jyi, an immediate application 
of Propositions 14.51 and 14.71 (applied to the boundary problem (14. 15)) rather than 
to (|4.ip) shows that 

HVU + \\D y ^ a \\ + \\D 2 y 4, a \\ < Ce 2 \\p\\ H i . 

To estimate Dy^j we will argue essentially as in the proof of Proposition ^. 71 One 
starts by proving interior estimates, which are obtained by feeding the test function 
ip := d 2 (x 2 d 2, ip) in the identity (|4.9p . As before, x(u) is a smooth function which 
vanishes in a neighborhood of <9B 2 and is identically equal to 1 in a disk of radius R 
and the script a (which is not summed) denotes any y direction. In order to get 
boundary estimates, one can start by noticing that the same argument we have used 
above to control derivatives with respect to a also works for dgip = y\ c^V' - yi diip- 
Indeed, the L 2 norm of the commutator 

g := A(d e tp) - 9e(AV0 

is also bounded by Ce||p|| as a consequence of Propositions I4.5H4.71 and besides 
dgip satisfies the boundary value problem 

A{d e $) = d e p + g inS^xD 2 , 

d v {dgip) = on S] x <9B 2 , / <9 e ^(a, y') dy' = . 

It then follows that the L 2 norm of the derivatives D^dg^p are bounded by Ce 2 ||/o||^fi . 
Hence one can now differentiate the equation = p with respect to r and iso- 
late dr"tp to show that the norm ||9j?-0||§i xAfJ in a neighborhood of the boundary is 
bounded by Ce 2 ||p||#i, as claimed. 

This proves the induction hypotheses for J+k = 1. For higher values of J+k, the 
idea is exactly the same. Derivatives with respect to a (or 6) are dealt with by taking 
derivatives directly in the equation and invoking the induction hypotheses. To 
estimates the highest derivative with respect to r, one combines interior estimates 
with test function 

^ = ^+ fe+1 ( X 2 ^+ fe+ ^) 
with the direct isolation of d^ +k+2 ip in the equation d^! +k (Aip — p) = 0. □ 

4.3. Pointwise estimates. Taking into account that D y ip = D y 4>, we can com- 
bine the results in the previous two subsections to obtain H k estimates for ip in 
which the derivatives with respect to the fast and slow variables are controlled in 
terms of different (and optimal) powers of e. To begin with, we can put together 
Propositions 14.31 and 14.61 and Theorem 14.81 to arrive at the following bounds: 

Theorem 4.9. The functions ip and ip satisfy the H k estimate 

< C k \\p\\ H u , U\\ Hk ^C k e 2 \\p\\ Hk . 

for any nonnegative integer k. 
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Remark 4.10. It is not hard to prove that these estimates are optimal with respect 
to the dependence on e. To have a rough idea of why this is true, the reader might 
want to consider the problem 
(4.16) 

i„ + = p in Sj x B 2 , d v ip = on §\ x9D 2 , / ipdady = 0, 

which can be understood as a simplified version of the problem (|4.ip . Here and in 
what follows, A y ip := d 2 ip+d 2 ip stands for the standard Laplacian in the variables y. 

For future reference, it is convenient to invoke the Sobolev embedding theorem 
and record the following estimate for the C k norm of the function ip: 

Theorem 4.11. The function ip satisfies the C k estimates 

UWc* < ChWpWhu , \\D y 1>\\ C k < C k e 2 \\p\\ Hk+3 

for any nonnegative integer k. 

Proof. It follows immediately from Theorem 14.91 upon noticing that D y ip = D y ip 
and using the Sobolev inequality H^Hc* *S C|Mlff fc + 2 - ^ 

Remark 4.12. It is clear that the same bounds we have proved for the function tp 
hold true if we assume that ip solves the model problem (I4.16|) instead of Aip = p. 
In particular, for future reference we record here that this function also satisfies the 
estimates 

(4.17) M\\ H >+' + ^ Dy ^ Hk < C k \\p\\ H u , \\D y iP\\ Hk ^C k e 2 \\p\\ H u. 

5. Harmonic fields in thin tubes 

In this section we use the estimates proved in Section U to compute the harmonic 
field in a thin tube up to terms that are suitably bounded for small e. 

We recall that a vector field h in the tube 71 is harmonic if it is divergence-free, 
irrotational, and tangent to the boundary. The vector space of all harmonic fields 
in the tube will be denoted by 

(5.1) H(T e ) := {he C*°°(0,M 3 ) :divh = 0, curl/i = 0, h-u = 0}. 

It is standard that the space H(T e ) is one-dimensional, as it is isomorphic to the 
first cohomology group of the tube with real coefficients. 

Let us consider the vector field in x D 2 defined by 

(5.2) h :=B- 2 (d a + Td e ). 

It can be readily checked that ho is irrotational (with respect to the metric (|3.3[) ) 
and tangent to the boundary. 

By the Hodge decomposition [IS], there is a function tp such that 

(5.3) h:=h + \7ifj 

is harmonic. Clearly this is the only harmonic vector field in x D 2 up to a 
multiplicative constant. By we are denoting the gradient of the function tp 
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with respect to the metric (|3.3|) . that is, 

ft; *\ V7 / ^a+T1pe ip r Aipg + e 2 r 2 rip a 

(5.4) V^ = — - 2 d a + - I d r + j^—y 2 Be. 

The fact that the field h is divergence-free implies that ip solves the Neumann 
boundary value problem 

Aip = p in §JxD 2 , d v ip = on S,} x <9B 2 , 

where 

(5.5) p := eB~ 3 r(rK sin 9 - n' cos 9) 

is minus the divergence of h (and, as such, satisfies J pdV = 0). We will also 
assume that J ip da dy = in order to determine tp uniquely. 

In the following two sections we will need some estimates for the harmonic field h 
(or, equivalently, for the function ip) that depend on the particular form of the 
source term p. These estimates are obtained in the following theorem, where we 
calculate, up to some controllable error, some derivatives of the function ip that we 
will use later on. To simplify the notation, we will write 0(e n ) for any function \ 
satisfying the bound ||x||c fe ^ C^e" for all k (and also for numbers whose absolute 
value is smaller than Ce n , but the meaning should be clear from the context). 

Theorem 5.1. Consider the functions 

(5.6) tpo :— (tk sin 9 — k'cos#) , 

8 

1 V 4 (Y 4 — 2r 2 ) 

(5.7) (ft := i ^(TK 2 sin26»~KK'cos26l), 

96 

which are obviously of order C(e 3 ) and C(e 4 ) 7 respectively. Then ip is related to 
these functions through the estimates 

i> = 0{e 2 ) , 
Dyip = Dycpo + G(e 4 ) , 
d e i> = dgtpo + d e <pi + 0(e 5 ) . 

Proof. Let p a := trijn sin 9 — n' cos9). It is easy to see that the function ip (a,y) 
is the only solution to the problem 

Ayipo = e 2 po in § ( ' x D 2 , d r ipo\ r =i = , / ipo dy = 

Jo 2 

for any value of the angle a. 

Consider now the function tpi := ip ~ ip Ql which obviously has zero normal de- 
rivative on x <9B 2 and has zero mean because so do tp and ipo . The Laplacian 
of ipi is given by 

Aipi=p- Aip Q = {p- p ) + {^p- - A(p, 

Using Eq. (|5.5[) . the first term in brackets can be easily shown to be 
p — po = 3e 2 nr 2 cos 0{tk sin 6 — k' cos 0) + C(e 3 ) , 
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while the second can be easily dealt with using the formula (|4. 2[) for the Laplacian: 
— ^ Aipo = - 1 cos 6 d r ipo I + C(e ) • 

Using the definition of ipo, this yields 

A^i - Pi + P2 + 0(e 3 ) , 

with 

13e 2 Kr 2 r . „„ , „„ n 
pi := [kt sin 20 — k cos 2u\ , 

e 2 (3 - Ur 2 )Kn' 
P2 := ~ 8 

of order 0(e 2 ), so the estimates we proved in Theorem 14 . 1 1 1 then ensure that 

xpi = 0(e 2 ) and D-ytpi = 0(e 4 ) . 
This shows that ip = 0{e 2 ) and D y ip = D y cp + C(e 4 ). 

To calculate D y ip up to C(e 5 ) we will consider two auxiliary functions (pi and f2- 
The function ipi (a, y) is the solution to the problem 

Aj,(y9i = e 2 pi in x D 2 , Qr^ilr^i = , / tpi dy = 0. 

The existence and uniqueness of the solution are standard given that the RHS 
satisfies 

pi dy = 

for all a. In fact, the solution can be computed in closed form using separation of 
variables, which readily yields the formula for ipi given in the statement (Eq. (|5.7p ). 

The function (f2 is the solution to the problem 

dlf 2 + = P2 in Sj x D 2 , a r ^ 2 |r=i = , J ip 2 dady = 0. 

Again, the existence and uniqueness of solutions is standard because the RHS sat- 
isfies the zero-mean condition 

P2 da dy = . 

(The reason why we are including derivatives with respect to a in the definition 
of ip2 but not in that of tpi is that the source term p2 has zero mean when averaged 
with respect to a and y but not when averaged in y only.) Obviously tp>2 does not 
depend on 0. 

We can now estimate "02 := ^i — Pi — P2 using the same argument we used 
with ipi. We start by noticing that, by construction, d r ip2 — when r = 1 and the 
integral J s± xD2 ip2 da dy is zero. The Laplacian of ip2 is 



A^2 = pi + P2 - Ay>i - Ap 2 + 0(e 3 ) 
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To pass to the third line we have used the expression of the Laplacian in polar 
coordinates (Eq. (|4.2p ) and of (p%, as well as the estimates for ip 2 that stem from 
Theorem 14.111 and Remark 14.121 Another application of Theorem 14.111 shows that 
D y i/j 2 =0(e 5 ). Since dgip2 = and 

1p = Ifio + if i + tp 2 + ^2 , 

the theorem follows. □ 



Our goal in this section is to construct a Beltrami field v in the thin tube T e , 
tangent to the boundary, whose harmonic part is a fixed harmonic field h. We will 
be particularly interested in the way the Beltrami field v is related to the harmonic 
field h as the parameter A tends to zero, since in the next section it will be crucial to 
have good estimates for this relation in order to compute some dynamical quantities 
of the field v. 

This section is divided in three parts. In the first subsection, we prove an exis- 
tence result for a boundary value problem for the curl operator in tubes in which 
we can prescribe the harmonic part of the solution (Corollary 16. 3|) . In the second 
subsection we will provide estimates for an auxiliary vector equation with constant 
coefficients (Proposition 16. 6p . which are used in the third subsection to prove the 
desired estimates, optimal in e, for the boundary value problem under consideration 
(Theorem [6J§. 

6.1. An existence result for the curl operator. Let us begin by making precise 
what we understand by the harmonic part of a vector field w in the thin tube T t 
that is tangent to the boundary. We will define its harmonic part to be the vector 



ll^llz^Cn) J% 

that is, its projection to the space of harmonic vector fields H(T e ), as introduced 
in Eq. (|5.1[) . In this subsection, we will denote by 



the usual H k norm of a vector field in the tube 7^Cl 3 and write L 2 (T € ) for H°(7l). 
Throughout, we will use the notation C € for positive constants, possibly not uni- 
formly bounded in the small parameter e, that may vary from line to line. 

In the following proposition we present the basic existence result that we will 
use to show the existence of Beltrami fields with prescribed harmonic part. The 
result is probably known to some experts but we have not found it in the literature. 
The proof relies on a duality argument for a suitable energy functional and the 
Fredholm alternative theorem. 

Proposition 6.1. Let f be an L 2 vector field in T e which is divergence-free. There 
is a countable subset of the real line without accumulation points such that, if the 
constant A does not belong to it, then the equation 

(6.1) curlw — Xw = / , divw = 



6. Beltrami fields with prescribed harmonic part 



field 
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has a unique H 1 solution w that is tangent to the boundary dTe and has zero 
harmonic part. 

Proof. Let us consider the Hilbert spaces of vector fields 

T := {F e L 2 (T e ,R 3 ) : divF = 0} , 

W := {w € i^Cr^R 3 ) : w is tangent to d%} , 

where the tangency condition is to be understood in terms of traces. It is well 
known |15) that, because of the tangency condition imposed on W, the H 1 norm is 
equivalent to the norm 

Hlvv : = II curl «; 1 1 2 L 2 {Te) + || divw\\ 2 L2(Te) + \\Pw\\ 2 L2{Te) 

in the sense that, for all w € W, 

IMIw 



a 



^ \\w\\h1(%) < Ce|M|w • 



Consider the scalar product E : W x W — > R associated to the norm || ■ ||w, 
given by 

E[w, u] := y ^ curl u; ■ curl it + div io div u + Vw ■ Vu^j dx . 

By the Riesz representation theorem, for any L 2 vector field F there is a unique 
wf € W such that 



(6.2) E[tOF,u] = y F-curludx 
for all tie W. 

We claim that, for any F E J 7 , the condition (|6.2|) is equivalent to demanding 
that wf be an H 1 solution to the equation 

(6.3) cutIwf = F , divu>F=0 

in the tube % that is tangent to the boundary and has zero harmonic part. One 
side of the implication is obvious, so we only need to prove that if wf satisfies the 
condition (|6.2p with F € J-, then it solves Eq. (|6.3[) . 

For this, let us take suitable choices of the field u in Eq. (|6.2|) . Letting u = Vwf 
be the harmonic part of wf, we obtain that Vwf is zero. To see that 6xvwf is zero 
too, it suffices to take u = V^, with ip being the solution to the boundary value 
problem 



Aip = div wf in %. , ^^19^=0, J ipdx = 0. 

Notice that, in order to show that the solution exists and Vip belongs to the 
space W, we need to use that wf is in H 1 and tangent to the boundary. 
Since Vwf and Aivwf are zero, the condition (|6.2p can then be written as 



(6.4) / cm\G-udx+ (G x v) ■ uda = 

for all u G W, where 

G := curl — F 
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and da denotes the induced surface measure on the boundary. Letting u vary over 
the space of smooth vector fields of compact support in 7^, we immediately infer 
that G must be irrotational. Now that we know the first term in (|6.4[) is always zero, 
we can take an arbitrary field u tangent to the boundary to derive that G x v must 
be zero, so that G is divergence-free, irrotational and orthogonal to the boundary. 
By the Hodge decomposition theorem, this ensures the field G is identically zero. 
Hence we have proved that wp is the only solution to the problem (|6.3j) . Notice 
that, as an immediate consequence of the equivalence of the H 1 norm to || • ||w on 
W, 

(6-5) \\w F \\ HH%) ^C € \\F\\ L2(%) 

for some constant that depends on e but not on F. 

Let us consider the operator K mapping a field F G T to its solution wp, which 
we regard as a linear map K : T — > T whose image lies in W. It is standard in 
view of the estimate (|6.5|) that the operator K is compact, so by the Fredholm 
alternative the equation 

(6.6) Kw — —w = f 

A 

has a unique solution for each / 6 T and any constant 1/A that does not belong to 
the spectrum of the adjoint operator K* . This spectrum is a bounded countable 
set that only accumulates at zero. 

Let w be the only solution to Eq. (|6.6|) with / = —Kf/X, for an arbitrary field / 
in J- . Since w — K(Xw + f) belongs to the image of K, from the above discussion 
it stems that the field w is tangent to the boundary, divergence-free and has zero 
harmonic part, and by the definition of K satisfies curl w — Xw + f. The proposition 
then follows. □ 

Remark 6.2. The proof of the proposition shows that the values of A for which 
there is not a unique solution to the Eq. (16.11) are given by the reciprocal of the 
eigenvalues of the compact operator K* that we defined in the proof. 

As a direct application of the previous proposition, we derive the following corol- 
lary, which gives the existence result for Beltrami fields with prescribed harmonic 
part that we will use later on. Of course, it is apparent that this corollary is totally 
different from the results that one has on compact manifolds without boundary, 
where the Beltrami equation has no solutions but for A belonging to a countable 
set and these solutions (i.e., the eigenfunctions of curl) are necessarily orthogonal 
to harmonic fields. 

Corollary 6.3. For any constant A not belonging to certain e-dependent countable 
set without accumulation points, there is a unique solution to the equation 

curl v = Xv , div v — 

in the tube 71 that is tangent to the boundary and whose harmonic part is Vv = h. 



Proof. It follows immediately from Proposition 16.11 by taking v = h + w, with w 
being the only solution to Eq. (|6 . 1 1) with / = Xh. □ 
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6.2. Estimates for an equation with constant coefficients. In this section 
we will prove estimates for a curl-type equation with constant coefficients on the 
domain xD 2 . This equation with constant coefficients is closely related to the 
Beltrami equation on the thin tube T e and will be used subsequently to estimate 
the difference between the Beltrami field v and its harmonic part h. As before, the 
natural coordinates on the domain §^ x D 2 will be denoted by (a,y), and we will 
consider polar coordinates (a, r, 9) where convenient. 

Vector fields on x B 2 are regarded as functions W : x D 2 — > M 3 , whose 
components are denoted by 

W = (W a ,W u W 2 ) . 

We will sometimes write W y := (Wi, Wb). On these vector fields we will consider 
the action of the differential operators 

Div W := d a W a + dxWx + d 2 W 2 , 
Curl, W := (W 2 - d 2 W x , - d a W 2 , d a Wx ~ ^ 



and we will say that a vector field W is tangent to the boundary if y ■ W y = (in the 
sense of traces) on the torus |y| = I. Obviously, Div and Curb are related to the 
standard divergence and curl operators through a rescaling, but for our purposes 
this form of the operators is more convenient. 

We will also consider the functional Q that maps each field W to the real number 



Q[W] := j W a dady. 



This is obviously related with the projection onto the space of "harmonic fields" 
associated to the operators Div and Curb, which is spanned by the constant field 
(1,0,0), so when Q[W] = we will simply say that the harmonic part of Q[v] is 
zero. All integrals are taken over x D 2 unless otherwise stated. To control the 
behavior of the vector fields, in addition to the usual H k norms we will consider 
the norms 

(6-7) \\Wf Hi :=||Wg|^+e 2 ||Wg|^. 

These norms should not be confused with the norm || • ||^i that we considered for 
scalar functions in Section [4] We will simply write || • || e for the corresponding 
L\ norm (k = 0), and reserve the notation || ■ || for the usual L 2 norm, where factors 
of e do not appear. 

In the following proposition we compute the lowest eigenvalue of a self-adjoint 
operator associated with Curb to show how the L 2 norm of a vector field W can 
be controlled using that of Curb W: 

Proposition 6.4. Let W be an H l vector field on E>\ x D 2 with Div W = that is 
tangent to the boundary and has zero harmonic part. Then 

HCurk^^k. 

e 

Proof. It is an easy consequence of [2 9) that Curb defines an unbounded self-adjoint 
operator on the Hilbert space 

% := {W 6 L 2 (S^xB 2 ,R 3 ) : Div W = 0, Q[W] = 0, W tangent to the boundary} , 
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endowed with the scalar product associated with the norm || • || e . The domain of 
this operator consists of the H 1 vector fields in H such that Curl e W is also in 
H. Hence, to prove the proposition it is enough to see that the eigenvalues of this 
operator, in absolute value, satisfy > C/e. 

Exploiting the symmetry of the equations (that is, rotation of the angles a and 
it is not hard to see that the eigenvalue equation 

Curl e W = fj,W , 

with W divergence-free, tangent to the boundary and with zero harmonic part, can 
be solved in closed form. Indeed, the symmetry ensures that the eigenfunctions can 
be chosen of the form 

W = c ma+im9 ( Vl {r) e a + v 2 (r) e r + v 3 (r) eg) , 

where n,m are integers, Vj(r) are functions of the radial variable and the unit 
vectors e a , e r , eg are defined in the obvious way: 

e a := (1, 0, 0) , e r :— (0, cos0, sm9) , eg := (0, — sin#, cos 9) . 

Using this expression, a tedious but straightforward computation shows that the 
smallest eigenvalues in absolute value are ±j[^ /e, where denotes the hrst pos- 
itive zero of the Bessel function J\. □ 

In the proof of the main result of this subsection we will need the following 
identity: 

Lemma 6.5. Let W be a vector field tangent to the boundary ofSj x D 2 . Then 

e 2 \\d a W y f + \\D y W y f + \\d a W a f + e- 2 \\D y W a \\ 2 = || Div W/|| 2 

+ || CiiT\ £ W\\ 2 e + [ \W v \ 2 dad9, 

where we are writing da d9 for the induced surface measure on x <9ID> 2 . 



Proof. One can easily check that 

||DivM/|| 2 + ||Curl e IU|| 2 = \\d a W a \\ 2 

+ £ (e 2 \\d a W t \\ 2 + + J2 \\d iWj f + S, 

V ' / i j=1 



where 



S := > J (diWi djWj - djWi diWj) da t 

i,jE{a,l,2} 
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We can now integrate by parts to write 

2 



S=-J2 I (VjWi diWj + yjW a d a Wj) da d9 

H\xdw- 

= J2 I W 3 d t {y 3 Wi)dad8 
V / WfdadO. 



i=i 



To pass to the second identity we have used that y ■ W y is zero on the boundary, 
so the second summand in the first integrand vanishes, and taken advantage of the 
fact that W is tangent to the boundary to integrate by parts a second time. □ 

In the following proposition, which is the main result in this subsection, we prove 
an estimate for the operators Div and Curl e that is "optimal" with respect to the 
small parameter e. In the proof of this version of the inequalities we will also derive 
another in which the RHS has one derivative less than the LHS, as is customary. 
However, the estimate we will need later on is the former, which is why it is the one 
that appears in the statement. It is not hard to check that the dependence on e of 
both estimates is sharp. 

Proposition 6.6. Let W be any vector field in §| x D 2 that is tangent to the 
boundary and satisfies the equation 

Curl £ W = F , Div W = p 

for a scalar function p and a vector field F. Then 

\\W a \\ Hk + \\D y W a \\ Hk sc C k (e\\F\\ H , + \\p\\ Hk + \Q[W]\) , 

\\d a W a \\ Hk ^C k (\\F\\ H , + \\p\\ Hk ), 

\\W v \\ Hk +e\\d a W v \\ Hk + \\D y W y \\ Hk < C k (\\F\\ Hk + \\p\\ Hk + \Q[W}\) 

for constants that depend on k but not on e. 

Proof. We can write the field W as the sum of three fields: 

The scalar function ip that appears in the second vector field is defined as the only 
solution to the Neumann boundary value problem 

(6.9) Bli, + ^ = p inSjxD 2 , a^lsjxao^O, / i>dady = 0, 

and it should be noticed that the third vector field (which corresponds to the 
harmonic part of W) is constant. As a consequence of these definitions and the 
properties of W, the field V is tangent to the boundary, has zero harmonic part 
(Q[V) = 0) and satisfies the equation 

(6.10) CwUV = F, DivV = 0. 
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The H k estimates stated in Remark 14.121 (Eg. (14.171) ). applied to the boundary 
problem (|6.9|) , provide suitable control of the second field that appears in Eq. ()6.8[) 
(that is, the "gradient" part) and its derivatives, as they show that 



IIMIh* + H^IU* + llDyda ^ Hk < c\\ P \\ Hk , 



\Dy1>\\ H > , \\Dyd a ll>\\ B H , H^VIU 



<C||p|| ff ». 
e 

The third field in Eq. (|6.8[) is trivial to control as it is constant. Therefore, to prove 
the proposition it is enough to derive suitable estimates for the field V . 

Hence, our goal is to show that the vector field V satisfies 

(6.11a) \\V a \\ Hk + \\D y V a \\ Hk < Ce||F|| ff; , 

(6.11b) \\d a V a \\ Hk ^C\\F\\ H u, 

(6.11C) \\Vy\\ Hk + e\\d a Vy\\ H H + \\DyVy\\ Hk C\\F\\ H H . 

For this we start by noticing that, when applied to Eq. (|6.10l) . the L 2 estimate 
proved in Proposition 16.41 yields ||F|| e < Ce||F|| e , or equivalently 

(6-12) ||Vy|<Ce||F|| e , ||V„||<C||F|| e . 

To prove the inequalities (|6.11l) . we start by using estimates for boundary traces 
and interpolation to control the term 

S ■= \Vy\ 2 = \\Vy\\ 2 L 2( S l xdI} 2\ 

< c\\v y \\ 2 L2(si)xH1/2{B2) < qi^iKii^n + iid^ii) 

that appears in Lemma 16.51 Together with the L 2 estimate (|6.12[) , we can then 
apply Lemma 16.51 to the field V to infer that 

\\d a V a \\^C\\F\\ e , \\D y V a \\^Ce\\F\\ e , 

\\d a Vy\\ < °\\F\\ e , \\DyVy\\ <C||F|| e . 

This proves the estimate (|6 . 1 1 1) for k = 0. 

The derivation of higher-order estimates from these inequalities is standard. To 
show the basic ideas, let us sketch the proof of the k = 1 estimates. Interior 
estimates are obtained by considering the identity shown in Lemma 16.51 for the 
vector field di{\V), where i = 1, 2 and % is a smooth function, compactly supported 
in E>] x D 2 and equal to one in the region \y\ < 1/2. Besides, it is easy to estimate 
derivatives of V with respect to a because the field d 3 a V is still tangent to the 
boundary, divergence- free, has zero harmonic part and satisfies the equation 

Curl e (^y) = d> a F. 

Analogous properties can be shown too for the globally defined field 

with de := yi 82V — J/2 d\ the angular derivative, so the same argument can be 
applied for this field. Since the second summand is obviously controlled by the 
previous estimates, this yields good H 1 estimates for dgV. To get estimates for the 
field V up to the boundary, it now suffices to control the radial derivative d r V in 
the region \y\ > 1/2, and this can be readily done by writing the equation (I6.10[) 
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in polar coordinates, taking its derivative with respect to r and isolating the terms 
having two radial derivatives. This process can be readily iterated to get the desired 
H k estimates. The details, which can be easily filled using these comments, are 
omitted. □ 

6.3. Estimates for Beltrami fields in thin tubes. In this subsection we es- 
tablish some estimates for Beltrami fields defined on the tube T e . The point of 
these estimates will be to control the difference between a Beltrami field v and its 
harmonic part h in terms of the Beltrami parameter A, while ensuring that the 
dependence on e of these estimates is optimal. As we did in Section [4] (but not in 
Subsection 16. we will identify T e with §J x D 2 through the coordinates (a, y). 

Given a vector field v in the tube, we will consider its components (v a , v±, v 2 ), 
defined through the relation 

v = v a d a + vi di + v 2 d 2 ■ 

We will use these coordinates to define the norms just as in Eq. (|6.7j) . of vector 
fields in T e , using the obvious formula: 

\H 2 H k := \\v a \\ 2 Hk +e 2 \\v y \\ 2 Hk . 

Here the Sobolcv norms in the RHS denote the usual scalar norms and we are using 
the notation v y :— (vi,V2). It should be noticed that these norms reflect to some 
extent the effect of the metric on the way vectors are measured; in particular, the 
L 2 norm defined by the above formula is obviously equivalent through constants 
that do not depend on e to the perhaps more appealing expression 

J g(v,v)dV, 

where g denotes the expression of the Euclidean metric in these coordinates and 
dV is the normalized volume (cf. Eq. (|3 . 1 1) and Eq. ([33])). 

We will begin with the following result, whose proof hinges on the analogous 
estimates proved for an equation with constant coefficients that we established in 
Proposition 16.61 

Proposition 6.7. Suppose that the vector field w is tangent to the boundary, has 
zero harmonic part (Vw = 0) and satisfies the equation 

(6.13) curlw = /, divw = 

in T e . Then 

\\w\\ H k sC Cfce||/|| H * . 

Proof. In the notation of Subsection 16. 2\ let us consider the vector field W : S\ x 
D 2 -> M 3 given by 

W := (Aw a + e 2 T(y 2 wi - yiw 2 ),wi + Ty 2 w on w 2 - Tyiw a ) , 

where we recall that A is the function defined in Eq. (|3.2D . It can be checked that W 
is given by the components (in the coordinates (a, y)) of the 1-form dual to w up 
to a rescaling of the y components by a factor of e~ 2 . It is worth noticing that, as 
the field w is tangent to the boundary of the tube T e , a simple computation using 
the definition of W shows that the field W is tangent to the boundary E>] x c© 2 , 
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that is, y ■ W y = on \y\ = 1. Moreover, the components of the field w on the tube 
can be recovered from those of W through the relations 



(6.14a) 
(6.14b) 
(6.14c) 



W a + e 2 T(y 1 W 2 -y 2 W 1 



Wi = W\ — 



w 2 = W 2 + 



where the function B was defined in Eq. (|3.6|) 
Using the notation 

B 



Si 



B 2 

ry 2 [W a + e 2 T(y 1 W 2 -y 2 W 1 )} 
B 2 

ryi[W a + e 2 r{ yi W 2 - yjWi)} 
B 2 



B 



Bo 



B 2 + e 2 r 2 y 2 ' ~ A ' B 2 + t 2 r 2 y\ 

we will also consider the field F : x ED 2 — > M 3 and the function p w : §J x D 2 
defined by 

F := 



Af a + e 2 T(y 2 f 1 - 7/1/2) fx + ry 2 f a f 2 - Tyif a 



PW '■= d a 



WJl 



B 
1 

B 



Bo 



wA 1 



1 

si 



do 



W 2 1 



1 

B~ 2 



Notice that, by the definition of the functions A and B, we trivially have 

(6.15) \\F\\ H} < C\\f\\ H{ , 

(6.16) \\pw\\h> < Ce(\\W\\ H * + \\d a W a \\ H * + \\D V W V \\ H >) . 

A straightforward computation shows that Eq. (|6.13p can be written in terms of 
the field W as 

(6.17) Curl e W = F , Div W = p w . 

To derive estimates for W, we start by analyzing Q[W]. To this end, we recall 
(cf. e.g. [2H]) that if da denotes the induced surface measure on the disk 

{a = a } C T e , 

the fact that the field w has zero harmonic part (Vw = 0) is equivalent to the 
assertion that 



0=/ g(w,v)da= Bw a \ 

J{a=a } JD 2 



dy 



for each angle ao- Since A — B is of order e, using the above equality one can 
estimate 



\Q[W}\ 



W a da dy 



(Aw a + e 2 T(y 2 wi - yiw 2 )) da dy 



J ((A - B)w a + e 2 r(y 2 w 1 - yiw 2 )) da dy 
< Ce|K|| +Ce 2 |K|| 



(6.18) < Ce|| W Q || + Ce 2 || W y || , 

where to derive the last inequality we have used the expression of the components 
of w in terms of those of W given in Eq. (|6.14p . 
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We are now ready to derive some estimates for the vector field W. Since W is 
tangent to E>j x cW 2 , we can apply Proposition 16.61 to Eq. (|6.17p to obtain 

(6.19a) \\W a \\ H n «S C(e\\F\\ H H + \\pw\\ H « + \Q[W]\) , 

(6.19b) \\d a W a \\ Hk < C(\\F\\ Hk + \\pwWhm) , 

(6.19c) \\W y \\ Hk + \\D y W y \\ Hk ^ C(\\F\\ H{ + \\pwWhx + \Q\W]\) ■ 

Going back to Eq. f|6 . , this yields 

\\ Pw \\ Hk ^ Ce(\\F\\ Hi + \\W a \\ H * + \\W y \\ Hk ) . 

Plugging this inequality and the bound (|6. 18|) for Q[W] in Eq. (|6 . 19[) . we immedi- 
ately get 

\\W a \\ H *^Ce(\\F\\ Hi + \\W a \\ H * + \\W y \\ Hk ), 
\\W y \\ Hk < C {\\F\\ H u + e\\W a \\ Hk + e\\W y \\ Hk ) . 
From these estimates one readily infers that, for small e, 
\\W a \\ Hk +e\\W y \\ Hk ^Ce\\F\\ Hk . 

In view of the formulas (|6.14l) . ^ CUVt^H^fc, so the proposition follows from 

the bound |6J5|) for the norm of F. □ 

We are now ready to prove the main result in this section, which estimates the 
difference between a Beltrami field in the tube and its harmonic part in terms of 
the Beltrami parameter A: 

Theorem 6.8. Let A be any nonzero real constant that is smaller in absolute value 
than some fixed positive constant A. For small enough e, the problem 

curlu = Xv 

has a unique solution in the tube l~ e that is tangent to the boundary and whose 
harmonic part Vv is the harmonic field h. Moreover, the difference between the 
field v and the harmonic field is bounded pointwise by 

\\v a - Mlc"=(S]xB> 2 ) + e \\ v y ~ h-y\\c k (Slx0 2 ) ^ Ck,A e \M ' 
where the constant only depends on k and A. 

Proof. Let us write w := v — h, so that the field w is tangent to the boundary, has 
zero harmonic part and satisfies the equation 

(6.20) curl w = \w + \h , divw = 0. 

We proved in Proposition 16.11 and Remark 16.21 that this equation has a unique 
solution if and only if 1/ A is not an eigenvalue of the operator K* introduced in 
the proof of the aforementioned proposition. As A is real, this is equivalent to 1/A 
being eigenvalue of the operator K , which means that there is a divergence- free 
L 2 field u such that 

A 

Suppose that 1/A is an eigenvalue of K. By the properties of K proved in 
Proposition 16.11 this means that the nonzero vector field Ku is tangent to the 
boundary, has zero harmonic part and satisfies 

cvx\{Ku) = XKu , div(Ku) = . 
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Applying the estimate proved in Proposition 16 . 71 to this equation, we get 

\\Ku\\ Ll ^Ce\X\\\Ku\\ Ll , 

which means that 

|A| > C/e 

for some positive constant that does not depend on e. 

Therefore the problem (|6.20p has a unique solution for all |A| ^ A, provided e is 
small enough. Moreover, when applied to Eq. (|6 . 20[) . Proposition 16 . 71 ensures that 

IMIjt* <C k e\X\(\\w\\ a} + \\h\\ Bi ), 

which yields 

\\w\\ H k < Cfce|A|||/i||fl-j 

provided that e is smaller than some constant of the form C/A. Since the norm 
of the harmonic field is bounded uniformly in e as 

due to Eq. (|5.3p and the estimates in Theorem 15. 11 we obtain 

|MU* < C k e\X\ . 

C s pointwise estimates for w are immediately obtained from this bound by taking 
k = s + 2 and using Sobolev embeddings. □ 

7. A KAM THEOREM FOR BELTRAMI FIELDS WITH SMALL A 

Let us consider the harmonic field h in the tube as introduced in Eq. (|5.3p . 
As before, we will assume that the thickness e is small. By Theorem 16.81 for any A 
smaller in absolute value than some fixed, e-independent constant, there is a unique 
solution v to the Beltrami equation 

curlu = Aw 

in % that is tangent to the boundary and whose harmonic part is h. We are 
interested in the case where the Beltrami constant A is suitably small. For simplicity 
of notation, throughout this section we will take A :— e 3 (although we could have 
taken any nonzero constant A = 0(e 3 )) and refer to the vector field v corresponding 
to this choice of A as the local Beltrami field. 

Our objective in this section is to study some fine dynamical properties of the 
local Beltrami field v. More precisely, we will show that for small values of e 
and "most" core curves 7, the boundary of the tube 71 is an invariant torus of 
the field v that is preserved (i.e., there is a small perturbation of dTe that is still 
invariant) under suitably small perturbations of v. For this, we will see that the 
key point is the analysis of the harmonic field h, which is close to the local Beltrami 
field v as a consequence of the estimates proved in Theorem 16.81 and the fact that 
the Beltrami parameter A = e 3 is small. At this point, it is worth emphasizing 
that, to some extent, the proofs of the dynamical properties that we study in this 
section ultimately depend on the estimates derived in Section 21 In particular, if 
the dependence on e of these estimates were worse, we would not be able to check 
the nondegeneracy conditions of the KAM theorem we prove in this section. 

This section is divided into four parts. In Subsection 17.11 we consider the tra- 
jectories of the local Beltrami field, after a rescaling of the field (which does not 
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alter their geometric structure), and calculate these trajectories perturbatively in 
a suitable range of time (Proposition l7.1[) . In Subsections 17.21 and 17.31 we use these 
perturbative expressions to compute the rotation number and normal torsion of 
the Poincare map of the local Beltrami held (Theorems 17.41 and 17. 8[) . which are 
the quantities that control the stability of the torus in the KAM theorem that we 
establish in Subsection 17.41 (Theorem I7.10[) . 

7.1. Trajectories of the local Beltrami field. In this subsection we aim to 
compute the trajectories of the local Beltrami field v perturbatively in the small 
parameter e. From Theorem l6.8l and the fact that the Beltrami parameter is A = e 3 , 
it follows that the local Beltrami field v is close to the harmonic field in the sense 
that 

(7.1) - ft. Q || C fc( S i xD 2) < C fe e 4 , \\v y - ^ y ||c fc (s«xD 2 ) < C k e 3 ■ 

Since the harmonic field can be written as h = ho+Vip (with ho and Vtp respectively 
given by (|5.2j) and (|5.4|1 ). from the estimates for ip proved in Theorem 15. II wc infer 
that 

v a =S- 2 (l + 7/; Q +r 1 /. e ) + 0(e 4 ) 
= l + 0(e). 

In particular, as v a does not vanish for small e, we can consider the analytic vector 
field 

(7.2) X := V 



l-'n 



We will use this vector field to study the geometric structure of the trajectories of 
the local Beltrami field, since both fields have the same unparametrized trajectories 
and the vector field X presents certain computational advantages, as we shall see 
in the next subsection. Before we go on, and identifying the tube T t with §J x B 2 
through the coordinates (a,y), let us note that the field X is well defined in a small 
neighborhood of Sj x D 2 because so is the local Beltrami field v. 

The trajectories of X are given by the parametrization (a(s),r(s),9(s)), where 
these functions satisfy the system of ODEs 

(7.3a) a = 1 , 

(7.3b) r = 9/i B ^ r — + 0{e 3 ) , 

(7.3c) . = T + ( £r) -^ + ^ 

1 + tpa + Tlpg 

These equations can be read off from the definition of the field X and its connection 
with the harmonic field (Eq. (|7.ip ) and the formulas (I5.3[) and (15.41) for the harmonic 
field and the gradient of ip. The point of the trajectory not only depends on the "flow 
parameter" s, but also on the initial conditions (ao,^Oj^o) at s — 0. Without loss 
of generality, in this section we will always take ao = 0, and make the dependence 
of the trajectory on (ro,&o) explicit by writing 

(a(s; ro, 9o),r(s; r ,0 Q ), 6{s; r ,6 )) 

when appropriate. 
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Throughout, it will be convenient to denote by 9 not only the angular coordinate 
in § , but also its lift to the real line. It should be noticed that the formulas we 
will give below are actually valid for the lifted coordinate too, which will be of use 
in the following subsection. 

In the following lemma we will compute the trajectory of the field X at time 
s G [0, £] up to a controllable error. We will assume that ro is bounded away from 
so that the trajectory cannot reach the coordinate singularity {r = 0} at any time 
s G [0, £]. This is convenient in view of the terms 1/r 2 that appear in the equations 
and is not a restriction for the applications that we have in mind, as we will be only 
concerned with initial conditions near the invariant torus {r = 1}. For simplicity, in 
this lemma we will abuse the notation and denote by 0(e 3 ) a quantity Q(ro, 6*o, s) 
that is uniformly bounded as 

\d* o d l eo Q(r ,6 ,s)\ <C kl ei 

for ro in any fixed compact set of the interval (0, 1] (which is the domain where 
polar coordinates define a diffeomorphism) , 9q G S 1 and s G [0, £], provided e is 
small enough. 

Proposition 7.1. Consider the solution to the system f|7.3|) with initial condition 
(0,7-0,00) an d r o > 0. At time s G [0,1], this solution is given by 

a(s;r o ,0 o ) = s, 
r(s;r ,9 ) = r + O(e), 

9{s; r ,9 ) = (O) (s) + e0 (1) (s) + eV 2 >(s) + 0(e 3 ) , 
where each quantity 9^' (s) = 6^' (s; ro, 9q) * s of order 0(1) and given by 

9^(s) :=9 Q + f r(s)ds, 



or 



1.2, 



. . 12-5r 2 f 3(r4 + 2r2-3)K(s) K (0) 
(s) := — -± / k(s) 2 t(s) ds + J cos 9 sin (o) (a) 



32 7 " v "' ' v " , ~ ' 64r 2 

„2\2 



_ (3-^KO) singocosg (0) (s) + (27 ' 5 ° r ° 8 Xf° )K(S) sin2^)( S ) 

+ (27 + l4-3h^(0f . n2 ^ 
384rQ 

Proof. Converting the ODEs (|7.3p into integral equations, it is clear that, for s G 
(7.4a) a(s;r ,9 ) = s, 

(7.4b) r(, nM . r „ + jT ?JI -^_ (0 + <,(«■) , 

(7.4c) »(„„,«„) = «„+ p + W-'^ + ^ ^ + c^). 

7o 1 + Wo. + i~Ve 
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In these equations all the functions under the integral signs are evaluated along the 
trajectories, i.e., at the point 

(7.5) a = s, r — r(s;r ,9 ) , 9 = 9(s; ro, #o) • 

Let us solve the equations perturbatively. We start by noticing that, as a con- 
sequence of the bounds for ijj derived in Thcorcni l5.ll (and its connection with the 
functions fo,tpi introduced in this theorem), the integrands can be expanded in e 
as 

r + {er)~ 2 A%l)0 + rip a _ ( Aip e 



= (r + ^ + r^)(l-^+0(e 3 )) 



1 + Ipa + T^e 

(7.6) = T + 4§ + c)u-'M. 



(7-7) m B * T - T = ^ + 0(e 2 ). 

e 2 (l + - 2 



Since ro > and s £ [0,£] (which allows us to control the effect of the denominator 
1/r 2 for small enough e), we immediately infer from Eqs. (|7.4bl) and (|7.4c[) that 

(7.8) r( S ;r Q ,9o)=r + 0(e), 9(s;r ,9 ) = 0<°> (s) + 0(e) . 

Of course, we define each function 9^ (s) as in the statement of the proposition. 

Now that we have the zeroth-order expression of the trajectories, we will next 
compute them up to second-order corrections. For convenience we will use the 
notation 

m . , r 3 -3r m . . 13(r 4 -2r 2 ) 

*o(r):= — , K l{ r):= Q( , 

for the dependence of tpo and tpi on r, respectively. We start with the analysis of 
the radial coordinate of the trajectories. Using again Eq. (|T. T[) and the zeroth-order 
estimates for the trajectories (|7.8[> . we derive that 



r(s;r ,9 ) = r + 



d r M^ro + Q(e),9^(s) + Q(e)) + 



^ 9 r M- S ,r Q 9^)) d - s + 0{e2) 



ds + 0{e 3 ) 



r + eTZ' (r ) \ [t(s)k(s) sm0 {o) (s) - k'(s) cos (o) («)] ds + C(e 2 ) 
d 



= r - eft (r ) / — (k(s) cos ^(5)) + 0(e 2 ) 
(7.9) =r + er( 1 )(s)+O(e 2 ), 



where 



-d) ( S ) := 3 ^ r °) [ K ( 5 ) cos 0(°> (s) - k(0) 



cos t) Q 



To pass to the second line we have used the mean value theorem and the obvious 
C k bound ifi — C(e 3 ), and to complete the calculation we have just plugged in 
the formulas for ip (Eq. (|5.6p ) and the definition of 9^(s). In a totally analogous 
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manner we can compute 9(s;ro,9o) up to C(e 2 ) 

Aip, 



= ^°)(s) + 



{erf 



0(e 3 



fryo(5,r ,fl(°>( a )) , 
(er ) 2 



(s,r o +O( e ),0(°>(s)+O(e)) 

ds + 0(e 2 



ds + 0{e 3 ) 



g(0)( a ) + fgoM f [ T (s) K (s) C o S 0(°)(s) + K '(s)sin^ o )(s)] 



(7.10) 



0(O)( s ) +£0 (l)( s ) +O ( £ 2 ) 



To complete the proof we need to calculate 9(s;r ,9o) up to 0(e 3 ). The proce- 
dure is as above but the computations are more tedious. We start by noticing that 
the term Aipg that appears in the integrand (|7.6[) can be written as 



(7.11) 



Aip e = dg(p + {dgipi - 2enr cos 9 dgtp ) +0{e 5 ) 



where we have used the estimates we established in Theorem 15.11 (recall that <pi 
was introduced in Eq. (|5.7|l ). Notice that the first summand is of order 0(e 3 ), while 
the term in brackets is of order C(e 4 ). 

We can now use the integral equation (|7.4c[) for the trajectories, the expan- 
sion (|7.6p for the integrand (together with (|7.11|l ). and the expression for the tra- 
jectories up to second order terms (Eqs. (|7.9[) and (|7.10p ) to get 



9(s;r ,9 )=e {0 Hs)+ f 

Jo € r 

dgipi — 2enrcos9 dgipo 



ds 



(s,r +erW +O(e 2 ).0(«) +c6(< 1 ) +C(e 2 )) 



ds + 0{e 3 ). 

For simplicity of notation, we omit the argument s when there is no risk of confusion. 
An elementary Taylor expansion and the mean value theorem show the first integral 
(let us call it h) is given by 



h = 



dg<po 



e9 (1 \s) + 



+4^V (i) (s)+%* (i) oo 



(s,r o ,0(°)) 



ds + 0{e 3 ) 



ds + 0(e 3 



(S,r ,e(0)) 



The second integral, which we call li , can be immediately simplified using the mean 
value theorem, finding that 



defi — 2enr cos9 dg(po 



ds + 0(e 3 ). 

(s,r o ,0(°>) 

The integrals I\ and 1% can be computed in closed form after replacing the 
functions tpo, tp± by their expressions, given in Eqs. (15.61) and (I5.7p . For example, 



h = Ji + .h + 0{e 3 ), 
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where 



Jl 



e 2 ^i(r ) 



e 2 ^iM 

2r o 
13 £ 2 (r§-2) 
96 



ds 

(s,r ,e(°>) 

I r(s>(s) 2 cos 2<9 (0) (s) + k(s)k'(s) sin 26» (0) (a)] 



^ K (s) 2 sm26^(s))ds 



[k(s) 2 sin 26l (0) (s)-K(0) 2 sin26»o] 



and 



J 2 := -2 



KCOSQ dg(fio 



2e 2 K (r ) 



ro 



ds 

s 

0(0) 

TKCOS 

0(0) 

+ k' sin 



e 2 ^ (r ) 



i-f K 2 sm20(°) 



rf.s 



e 2 (3 



16 



1 + k(s) 2 sin 20 (o) (s) - k(0) 2 sin 26 



The other terms can be dealt with using analogous arguments, arriving at the 
formula for 8^ (s) that appears in the statement. □ 

To conclude this subsection, we will show that the trajectories of the field X on 
the invariant torus r = 1 satisfy certain functional equation up to some controllable 
errors. The reason why we need to consider this way of describing trajectories on 
the invariant torus is that, in order to compute the rotation number later on, we 
will need to understand the trajectories on the torus for arbitrarily large times. 
The expression for the trajectories we obtained in Proposition 17. II is not well suited 
for this purpose, while the functional equation below turns out to be much more 
convenient. In order to describe the errors that appear in the functional equation, in 
the following proposition we will use the notation s 0(e n ) for any quantity Q(9q, s) 
that is bounded as 



KQ(6o,s)\ <C,-(l + | S |)e" 



and \d a 8LQ(6o,s)\ < Cj€ n 



for nonnegative integer j (we could have considered higher derivatives with respect 
to s too, but we will not need this feature). 

Proposition 7.2. Consider the trajectories of the system of ODEs (|7.3|) with ini- 
tial condition (cto,ro,6o) = (0,1, 9q)- The function 9(s) = #(s;l,#o) satisfies the 
approximate functional equation 



+ / r(s)ds-- k(s) sin 6»(s) - k(0) sin O ] +sO(e 2 ). 



Proof. The starting point is the differential equations for the trajectories (|7.3[) with 
initial radius ro = 1. It is obvious that the radial component of the trajectory is 

r(s;l,e ) = 1, 
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which simply shows that the set {r = 1} is an invariant torus. Therefore, Eq. (|7.3c|) 
for 9(s) becomes 



(7.12) 



T + € 2 Alpg + tV>o 



2 -Aip 6 - T 2 il> e 



+ 0(e 3 ) , 



1 + Ipa + Tlpe 1 + 1p a + T1p0 

where the functions in the RHS are evaluated along the trajectories (s,l,9(s)). 
Expanding the fraction that appears in the second identity using the estimates 
for ip proved in Theorem 15. II we therefore arrive at 

(7.13) 



■ 5) = T(s) + ^(^ +0(e2) 



Notice that the fraction is of order 0(e). Converting (|7.12[) into an integral equation, 
an immediate consequence of this estimate is that 

j s d fl y> (a, 1.0 (a)) 

Jo f2 



9{s) = 0(°>(* 



-ds + sO(e 2 



This integral can be evaluated, modulo sO(e 2 ), using the formula ()7.13|) for the 
derivative 9, the expressions of tpo, and integration by parts: 

Jo « 2 4 . 

e 



HT COS f 



k sin 



9] ds 



k(9 + £>(e)) cos 9 + «' sin* 



ds 



k sin 9 )ds + sOie 2 



d 

4J ds 

= ~[k(s) sin 9(s) - k(0) sin 9 ] +sO{e 2 ). 
Here, of course, all the integrands are evaluated at the point (s, l,9(s)). 



□ 



7.2. Rotation number of the Poincare map of the local Beltrami field. We 

will denote by <f> s the time-s flow of the field X, which maps each point (aoj r 0i $o) 
to the trajectory of the ODEs (|7.3|) at time s that has the latter values as initial 
conditions. Since the field X, introduced in Eq. (|7.2j) . is tangent to the boundary of 
the domain Si x B 2 , it is standard that the flow 4> s is a well defined diffeomorphism 



of § 



J 2 for all values of s. 



Let us now consider the Poincare map of the field X, which is the tool we will use 
to analyze the dynamical properties of the flow (and which coincides with that of 
the local Beltrami field v). For this, we start by considering the section {a — ao}, 
which is clearly transverse to the vector field X. The Poincare map of this section, 
n Qo : D 2 — > D 2 , sends each point (ro,0o) <= ^) 2 to the first point at which the 
trajectory <fi s (ao,ro,9o) intersects the section {a — ao} (with s > 0). The reason 
why we are considering the field X is that it is isochronous in the sense that this 
first return point is given by the time--^ flow of X, that is, 



(7.14) 



n ao (r ,6>o) = 0i(a o ,r o ,9 o ) . 



We will omit the subscript when ao — 0, and use Cartesian coordinates y in the 
disk when convenient. 

One should notice that, since we are assuming that the curve 7 is analytic, the 
boundary of the tube % is also an analytic surface, so it is standard [53] that 
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the field v is analytic in a neighborhood of the closure T e - This ensures that the 
Poincare map is also a well-defined analytic map in a neighborhood of the closed 
disk W. 

In the following proposition we will show that the Poincare map of the harmonic 
field preserves a measure on the disk. For later convenience, we will state this result 
in terms of the associated 2-form rather than the measure: 

Proposition 7.3. The Poincare map II preserves the positive measure on the disk 
corresponding to the 2-form 

A := G A (r, 6)rdrAd6 

on the disk D 2 , with 

(7.15) G A (r,6) :=Bv a \ a=0 = 1 + e«(0) r cos 9 + 0(e 2 ) . 

Proof. That the function G\(y) := Bv a \ a= o has indeed the form given by the RHS 
of (|7.15[) is an immediate consequence of the estimates for the function ip proved 
in Theorem l5.ll and Eq. (|7.ip . Given a Borel set 6cD 2 and a small positive 6, let 
us denote by 



:= / A 
Jb 



its area and let 

B 5 := {-6,5) x B 

be a small thickening of the set {0} x B in the closed domain §J x B 2 . 

Since the divergence of v is zero, from the definition of X it stems that its flow 
preserves the volume 

dV := v a dV . 
Clearly the V- volume of the set Bs is 

V(B 5 ):= [ dV = [ [ Bv a dyda 
JBs J-8 Jb 

(7.16) =26(1 + 0(6)) f Bv a \ a=0 dy = 26fi(B) + O(6 2 ). 

JB 

Let us now observe that the image of the set Bs under the time-£ flow fa is given 

by 

(7.17) MBb)= IJ W xn »( B )' 

~8<a<S 

By the continuous dependence of the flow on the initial conditions and Eq. (|7.14l) , 
the Poincare maps corresponding to different values of the angle a satisfy 



l|n Q - ii||c°(b 2 ) ^ C8 
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for \a\ < 5, so we can use the decomposition (|7.17|) to show that the V- volume of 

4>i(Bs) is 

V(<t>t(B s )) =11 Bv a dyda 
J -s Jn a (B) 

c6 



L(jL*- | -* +0 " ) )*' 



KB) 

(7.18) = 26^(U(B)) + 0(S 2 ). 

Equating the F-volumes of Bs and (j)g(Bs), given by Eqs. (|7. 16|) and (|7. 18[) . and 
considering small values of S we then obtain that 

l*(B) = ^U(B)) , 

as claimed. □ 

Since the local Beltrami field v is tangent to the boundary of the domain SjxD 2 , 
the image of dD 2 under the Poincare map II is also contained in c© 2 . Hence, the 
restriction of II to <90 2 defines an analytic diffeomorphism of the circle, which will 
be denoted by 

n| aD 2 : <9B 2 -> do 2 . 

Using the coordinate 9 to identify the circle dti 2 with M/27rZ, the latter circle 
diffeomorphism can be naturally lifted to a diffeomorphism of the real line that we 
will denote by II : 1 -> 1. As is well known, a basic tool in the study of circle 
diffcomorphisms is the rotation number (or frequency) of the map, which is defined 
as 

(7.19) Wn:= i im Mj5. 

n—^oo n 

Here ft" denotes the n th iterate of II and Oq is any real number. Since the Poincare 
map of a flow is homotopic to the identity, it is standard that the above limit exists 
and is independent of the choice of 9q [28] . 

We shall next compute the rotation number of the circle diffeomorphism II|aD 2 
using the functional equation satisfied (up to controllable errors) by the trajectories 
of the field X on the invariant torus. The reason is that, in order to compute the 
rotation number to order C(e 2 ), we need to iterate the Poincare map an arbitrarily 
large number of times, which requires fine control of the growth of the errors for 
large times. 

The following theorem asserts that the rotation number is given by the total 
torsion not only modulo 0(e), as can be shown without relying on the functional 
equation, but also modulo C(e 2 ). The fact that the 0(e) correction is zero will be 
important later on. 

Theorem 7.4. The rotation number of the circle diffeomorphism H\go2 is 

u n = I ' T(a)da + 0(e 2 ). 
Jq 

Proof. By the definition of the flow, Eq. (|7. 19[) simply asserts that 

9(n£)-9 
luu := lim , 
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where 8(n£) = 9(n£; 1,0q) denotes the angular component of the trajectory solving 
the system (17. 3[) with initial condition (0,1, 0o), evaluated at time n£. Since the 
curvature n(a) and torsion r(a) are ^-periodic, Proposition 17.21 then ensures that 

rds-^ [sin 0(nl) - sin go] + nlO(e 2 ) 
wn = hm — d = / ras + 0[e ) . 

n-s-oo n J 

□ 

7.3. The nondegeneracy condition for the Poincare map. In this subsection 
we will compute a quantity associated with the Poincare map II (sometimes called 
the normal torsion of the map) that was introduced to analyze the stability of 
individual invariant tori of symplectic diffeomorphisms [HI \E\ ■ As we shall see, the 
assumption that the normal torsion is nonzero plays a role that is analogous to the 
twist condition in the classical theorem by Arnold and Moser on perturbations of 
intcgrable symplectic maps. As the name can be misleading, it is worth emphasizing 
that, in principle, the normal torsion has nothing to do with the torsion of a curve. 

Let us begin by introducing some notation. We will consider a domain T> in the 
plane that contains the closed unit disk D 2 and a map II : T> — > M 2 . (Eventually, 
we will be interested in taking as II the Poincare map II introduced in the previous 
subsection.) A closed curve r C T> is invariant if its image II(L) is contained in 
r. If r is an invariant curve of II, one says that II|r is conjugate to a rotation of 
frequency uj through the diffeomorphism : S 1 — > L if 

e _1 o n| r o e(i?) = $ + uj, 

for all § in S 1 . When L = <9D 2 , we will abuse the notation and also denote by 
the diffeomorphism S 1 — > S 1 corresponding to the angular component of the above 
diffeomorphism S 1 — > dU 2 . (Therefore, in the case of dH) 2 the above diffeomorphism 
will read as y = (cos 0($), sin 0(i?)) in Cartesian coordinates and (r, 0) = (1, ©(#)) 
in polar coordinates.) From the context it will be clear which interpretation of 
must be considered in each case. 

Definition 7.5. Let II : D 2 — ► D 2 be a diffeomorphism of the disk that preserves 
the measure defined by the 2-form G(r, 9)rdr A d6. We will denote the radial 
and angular components of II by (LI r , lie), respectively. Assume that n|g D 2 is a 
diffeomorphism of the circle <9B 2 — > 3D 2 that is conjugate to a rotation of frequency 
io through a diffeomorphism 0, which we regard here as a map S 1 — > S 1 . The normal 
torsion of the map LI on the invariant circle <9B 2 is the real number 



N - = r - wi,0 W) 

n ' ' e'OJ + wje'O?) G(i, ©(<?)) 



The reason why we consider the above quantity is that it appears in a nondegen- 
eracy condition of a theorem by de la Llave et al. [8]. In fact, the result [8j Theorem 
1] is much more general, and we will only need a concrete application that we state 
next in a form that is particularly well suited for our purposes. The normal torsion 
was also considered for the same purpose by Herman in |18j when G(r, 9) = 1 and 
0(i?) = Before stating the theorem, let us recall that a number to is Diophantine 
if there exist a positive constant C and v > 1 such that 

(7.20) " C 



uj p 
2^ k 
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for any integers p, k with k 1. 

Theorem 7.6 (De la Llave et al. [8]). Consider a small neighborhood T> of the closed 
unit disk D 2 in R 2 . Take an analytic map II : T> — > R 2 that is a diffeomorphism 
onto its image preserving the measure G(r, 9) r dr d9, with G analytic in T> . Suppose 
that the following two conditions hold: 

(i) The circle <9B 2 is invariant, and Elian 2 is conjugate through an analytic 
diffeomorphism : S 1 — > SO 2 to a rotation whose frequency lo satisfies a 
Diophantine condition. 

(ii) The normal torsion J\f^ of the map II on the invariant circle dU> 2 is 
nonzero. 

Then for each 5 > and positive integer m there are 5' > and an integer k such 
that, if an analytic map II : D — > R 2 preserving the same measure G{r, 6) r dr d9 
satisfies 

\\n-n\\ ck{v) <5' , 

one can transform the circle 9D 2 by a diffeomorphism ^> of R 2 so that ^(cTO 2 ) is 
an invariant curve for the map II. Moreover, the map H\^fdjp) is also conjugate to 
a rotation of frequency uj and the difference ^ — id can be assumed to be supported 
in a small neighborhood of dU> 2 and such that 

||*-id|| c ™ 

Proof. The statement is simply a rewording of J8J Theorem 1] , in the particular case 
of planar maps and omitting some quantitative estimates that will not be needed 
in the rest of the paper. The only point that requires more elaboration is to check 
what the degeneracy condition looks like in the situation we are considering in this 
section. 

For the benefit of the reader, let us give some details about how the statement is 
derived from [HI Theorem 1], borrowing some notation from this reference without 
further mention. The map II is obviously symplectic with respect to the analytic 
2-form A := G(r,9) r dr A d6. This 2-form is obviously exact, as A = d/3 with the 
smooth 1-form on T> 



P : =[J G A (f,6)fdfj de. 

Moreover, the 1-forms IT/3 — j3 and U* (3 — f3 are exact because 

(fl*/3-/3)= / (dft* f3 - d/3) = /(ft*A-A) = 0, 



and analogously for II. The maps II, II and G being analytic, it is clear that, by 
taking a smaller neighborhood of D 2 if necessary, one can holomorphically extend 
these maps to a small neighborhood in C 2 of T> C R 2 C C 2 . 

As before, let us now regard as a diffeomorphism of S 1 . Consider the embed- 
ding K : S 1 — > T> given by 

K(-&) := (cos 8(0), sin 0(0)) , 

which is analytic by hypothesis. By periodicity, the embedding K can be holomor- 
phically extended to a bounded function on a thin complex strip {$ G C : | Im0| < 
c}. 
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Therefore, the only hypothesis of [H Theorem 1] that is not immediate is the 
non-degeneracy condition. Let us take Cartesian components (yi, j/2) in I? and call 
(III, II2) the Cartesian components of the map LT. In the aforementioned reference, 
the condition is that the average of the function 

S{$) := P(0 + u) ■ [DU(<&) J(0) -1 P(0)] 

be nonzero, where the dot denotes the Euclidean scalar product. Here P($) := 
&('&)- 1 (-8m@('d),cose('d)), 

(Ufi) y (i?) :=^(JfW) 

is the Jacobian matrix of LT and 

J(tf) := G(JjT(0)) (J . 

Expressing the Cartesian components in polar coordinates, 
III = LT r cos lie , II 2 = II r sin tig , 
we immediately obtain that 

_ d r tl r sin(f[e - 6(t? + u)) + gr^gg cos(n 9 - 9(i? + op) 

where the functions whose argument has not been specified are evaluated at K^). 
Since Tl r (K( , d)) = 1 because the circle <9B 2 is invariant and, by hypothesis, II|aD 2 
is conjugate to the rotation of frequency ui through the diffeomorphism (i.e., 

for all 1?), we get 

rr/ Q\ = 1 9 r n fl (^(^)) 

1 J 0'(i? + u;)0'(tf) G(ir(0)) ' 

In view of the way we defined the normal torsion (cf. Definition I7.5p . the statement 
then follows immediately from [8j Theorem 1] after realizing that 

□ 



In order to calculate the normal torsion of the Poincare map LT, let us begin by 
computing the diffeomorphism that conjugates II|g]j 2 to a rotation when its rotation 
number satisfies a Diophantine condition: 

Proposition 7.7. Suppose that the rotation number wn of the Poincare map of 
the local Beltrami field is Diophantine. Then the circle diffeomorphism n|an 2 is 
conjugate to a rotation of frequency wn through an analytic diffeomorphism : 
S 1 ->■ 5IB> 2 that satisfies 



9(0) = $ - IJp sm tf + o(e 2 ) 
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Proof. As the rotation number wn satisfies a Diophantine condition, the map n|aij2 
is conjugate to a rotation of frequency wn through an analytic diffeomorphism |28[ 
Theorem 1.3]. This diffeomorphism can be understood as a change of coordinates 

e a = e(i?). 

Let us now compute the diffeomorphism 0. Writing 9o = 0(#), the fact that 
Il|do 2 is conjugate to a rotation of frequency u>n through means that the trajec- 
tory 9(1; 1, 0(i?)) at time i corresponds to 0(i? + wn), for any choice of i?. Using 
the equation for the trajectory proved in Proposition [721 at time £, this means that 
must satisfy the equation 

(7.21) 0(t?+wn) = 0(i?)+wn + ^^[sin0(i?)-sin0(i? + w n )] + 0(e 2 ) . 

Here we have used the expression for ujn proved in Theorem I7.4I (which allows us 
to replace f r ds = wn + 0(e 2 )) and that the function k is ^-periodic. 

To analyze this equation, let us write the 0(e 2 ) term in the RHS of (|7.21[) 
as R(d). Let us identify <9ID> 2 with S 1 = R/2-7rZ through the angular coordinate 
#o, so that is regarded as a diffeomorphism S 1 — > § . With a slight abuse of 
notation, let us still denote by its lift K — > R. It is standard that this lift can be 
written as 

(7.22) e(0)=0 + -H"(0), 
where H(-d) is an analytic 27r-periodic function. 

Let us consider the 27r-periodic function 

(7.23) :=If(tf) + ^^ sin ©(<?). 
Eq. (|7.2ip then reads as 

(7.24) F(i? + wn) = F(i?) + , 

with R = 0(e 2 ). Consider the Fourier series of the functions F and R: 

oo oo 

f(<?)= J2 p * eikS > R w= E 

k— — oo k— — oo 

Eq. (|7.24l) then asserts that for any nonzero integer k the Fourier coefficients of F 
and i? are related through the identity 

(7.25) F k ,h 



We can obviously take Fo = 0; moreover, Rq = because it is a necessary condition 
for the existence of the diffeomorphism 0. 

Since wn satisfies the Diophantine condition (|7.20|) , for large integer values of k 
we have the elementary inequality 

| e ifcc n _!| > C\k\-", 

so that from Eq. (|7.25[) the H m norm of F can be estimated by 

\\f\\h~ <c( £ (i + fc 2 ) m |fc| 2 Hi? fc | 2 ) 2 < a||ii||ff«+. < c m e 2 
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for any nonnegative integer m. To derive the last inequality, which shows that 
F = 0(e 2 ), we have used that R = £>(e 2 ). In view of Eqs. ([732]) and (f7T23|) . this 
ensures that 

&(§)=§- sin 0(0) + 0(e 2 ) . 

In turn, this readily leads to the expression for the diffeomorphism provided in 
the statement. □ 

We are ready to provide a closed formula for the normal torsion of the Poincare 
map, up to terms of order C(e 3 ). The leading term only depends on the geometry 
of the curve (through its curvature and torsion) and, as is to be expected, not on 
the section of S} x D 2 we used to define the Poincare map: 

Theorem 7.8. Suppose that the rotation number lou satisfies a Diophantine con- 
dition. Then the normal torsion of the Poincare map of the local Beltrami field v 
on the invariant circle 3B 2 is 

Ml = -^r / K(a) 2 T(a)da + 0(e 3 ). 
8 Jo 

Proof. By definition, the normal torsion is 

si- . = r 1 Wi.ew) 
n " Jo e>(# + m)&(d) G A (i,ew) ' 

where Ug (r, 9) is the angular component of the Poincare map, (i?) is the diffeomor- 
phism defined in Proposition 17.71 (considered as a map S 1 — > S 1 ) and the function 
G A (r, 9) is given by Eq. ([7TT5]) . 

We have already computed all the terms we need to evaluate the integrand up 
to an C(e 3 ) error. Indeed, from Propositions 17.11 17.31 and 17.71 and Theorem 17.41 it 
stems that, setting T := J r(a) da, 

oj n =T + 0(e 2 ), 

&(0 + wn) = 1 - cos(tf + T) + 0(e 2 ) , 

G A (1, 0(i?)) = 1 + e/e(0) cos-d + 0(e 2 ) 

«n« ^_de( 1] f,i * ^(o) a , 



Or 



■(1, 9 (tf)) = e^- (l; 1, * - ™ sin,?) + 1, *) + <V) , 



the functions 9^'(s;ro,9o) being in Proposition 17. II Plugging these expressions in 
the integral for the normal torsion and using trigonometric identities, one arrives 
at the expression 

^ r <M(0), . .. . a s he 2 ^ 



Mi 



'sin(i? + T) -sini?) / «;(a) 2 r(a) da 

>ir * e2 smTcos(2$ + T) + 0(e 3 ) 
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di) 



5ttc ' 



\ K 2 rda + 0(e 3 ), 
Jo 
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as claimed. □ 

7.4. A KAM theorem for generic tubes. In this subsection we will use the pre- 
vious results to prove a theorem on the preservation of invariant tori for divergence- 
free vector fields that are close to the local Beltrami field v for small enough e. As 
we shall see, the hypotheses of the KAM theorem will hold true as long as the core 
curve 7 of the tube satisfies certain generic geometric conditions. Details on the 
validity of these conditions are given below. 

Since the local Beltrami field v is analytic in a neighborhood of the closure of 
the tube, which we identify with §^ x D 2 via the coordinates (a, y), we can assume 
that v is defined in some domain x T>, with V a neighborhood of the closed 
unit disk in the plane. To measure the smallness of a field, we will use the norm 
ll' u llc fc (s 1 xD)i which we define in terms of its components in the coordinates (a, y) 
as 

\\ u \\c k (S]xV) '■= \\ u a\\c k ($ 1 e xV) + \\ u v\\c($\xT>) ■ 

We will sometimes find it convenient to refer to a domain bounded by an invariant 
torus as an invariant tube of the field. We say that a field u is orbitally conjugate 
to a rotation of frequency uj on an invariant torus S if there are global coordinates 
(a, 9) : £ —> S 1 x S 1 in which the vector field is linear with frequency a; up to a 
multiplicative factor, that is, 

u| B = F(a,0) (dz + udg), 

where F : S 1 x 8 1 — ► K is a nonvanishing function. 

In the following lemma we will show that, for a generic core curve 7, the rotation 
number of the Poincare map of the local Beltrami field is Diophantine and its normal 
torsion is nonzero. To make precise what we understand by "generic", we will say 
that certain property holds for a C m -dense set of closed analytic curves if, given 
any closed analytic curve 7 in M 3 , one can deform it by a diffcomorphism $ of K 3 , 
with ||$ — id|| Cm ( R 3) as small as one wishes, so that the curve $(7) has the desired 
property. Notice that this does not imply that the property holds for an open set 
of curves. 

Lemma 7.9. Let m be any positive integer. The set of closed analytic curves for 
which the Poincare map of the local Beltrami field v has a Diophantine rotation 
number Uji and nonzero normal torsion Afn on the invariant circle <91D> 2 is C™ - 
dense. 

Proof. The result is not hard to prove using the expressions for the rotation number 
and the normal torsion derived in Theorems 17.41 and 17.81 

r e 

un= / rda + 0(e 2 ) , 
Jo 

-^n = -—3— f K 2 Tda + (D(e 3 ). 
8 Jo 

A way of making things precise is the following. We will consider deformations 
of the curve 7, labeled by a parameter 8. More concretely, let us denote by 
(T(a),N(a),B(a)) the Frenet trihedron of the curve 7 at the point 7(a) (one 
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should not mistake the binormal vector B{a) for the function B that we intro- 
duced in Eq. (|3.6p . which will not be used in this proof). Let F (a) be an analytic 
^-periodic function and consider the family of curves r(a,<5) in M 3 given by 

T(a,5) := j{a) + 5 F{a)B{a) . 

If 6 is close to zero, 75 = T(-,(5) is a closed analytic curve. Notice that, for 6^0, 
a is no longer an arc-length parametrization of 75 but, due to the properties of the 
binormal field, 

|r'| = i + o(s 2 ). 

Here and in what follows, we denote by a prime the derivatives with respect to a. 
In particular, the length of 75 is £ + 0(S 2 ). We will label the geometric quantities 
associated with the curve 7,5 with a subscript S (e.g., k§ and t$ for its curvature 
and torsion). 

The results we have presented in this section carry over immediately when one 
does not only consider the tube T e associated with the curve 7, but the family of 
tubes Te("js)- The dependence of the various quantities on the small parameter S 
is smooth and can be controlled easily. A tedious but straightforward computation 
using the well-known formulas for the curvature and torsion of a parametrized curve 
shows that 



2Sk(2F't + Ft') + <D(S 2 ) 



TS = T - S 



kF' + — 



d ( F" — t F 



da\ k 



0{8 2 



Since the length of 75 differs from £ by an 0(S 2 ) term, one readily finds that the 
rotation number of the Poincare map associated with the harmonic field of the tube 
7^(7<5) is 



w ni = / T S (a) da + 0(e 2 + S 2 ) 
Jo 

= uj n + S [ n'(a) F(a) da + 0(e 2 + S 2 ) 







It is not difficult to see that the derivative of the above term C(e 2 + S 2 ) with 
respect to the parameter 5 is of order 0(e 2 + 6), so it follows that the derivative of 
the rotation number with respect to 5 at is 



(7.26) duJU 



,IS 



s=o Jo 



K,'(a)F(a) da + 0(e 



2\ 



Similarly, the variation of the normal torsion with respect to d can be shown to be 

57T6 2 



(7.27) dNn > 



dS 



— / (2/e'" + 3kV -6ktt' -6r 2 K ')Fda + 0(e 3 ). 
8 Jo 



6=0 

Perturbing the curve 7 a little if necessary to ensure that the functions k' and 

2k'" + 3k 2 k! - 6ktt' - 6rV 

are not identically zero, we deduce that the function F can be chosen so that 
the derivatives (|7.26p and (|7.27l) are nonzero. Since Diophantine numbers have 
full Lebesgue measure, this immediately implies that one can choose an arbitrarily 
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small 6 such that the rotation number ujyi s satisfies a Diophantine condition and 
the normal torsion 7Vn s is nonzero. □ 

We can now show that, for a generic core curve 7 and small enough e, a suitably 
small perturbation of the local Beltrami field still has an invariant torus that is 
close to the original one: 

Theorem 7.10. For any positive integer m there is a C m -dense set of closed 
analytic curves 7 with the following KAM-type property: for any S > 0, there is 
another positive integer k and some 5' > such that any analytic divergence-free 
vector field u whose difference with the local Beltrami field v of the tube T t =T t (7) 
is bounded by 

(7-28) ||u-«||c*(s}xx>) 

possesses an invariant tube. Furthermore, one can find a diffeomorphism ^> ofM 3 , 
with || "J — id||(7m( R 3) < 5 and $ — id supported in a small neighborhood of the 
torus dT e , such that ^(71) is an invariant tube of u and u is orbitally conjugate on 
the invariant torus d^f(7l) to a Diophantine rotation. 

Proof. By Lemma 17. 9[ we can deform the curve that lies at the core of the tube 
% by a diffeomorphism of R 3 , arbitrarily close to the identity in the C m norm, so 
that, if we consider the local Beltrami field (which we still denote by v) associated 
to the deformed tube, its rotation number wn satisfies a Diophantine condition and 
its normal torsion An is nonzero. As before, we will identify this deformed tube 
with the domain §) x D 2 through adapted coordinates (a, y). 

Consider the Poincare map n of the local Beltrami field v, which can be safely 
considered as a diffeomorphism n from T> onto its image, with T> a neighborhood 
of the closed unit disk D 2 . The Poincare map of the field u, also defined on the 
section {a = 0}, is another diffeomorphism of T> onto its image that we will denote 
by EL The vector fields u and v being close by (|7.28j) . it is apparent that 

(7.29) \\U-il\\ CHv) <CS' 

as the Poincare map is simply obtained by integrating the associated vector field 
along a trajectory between two consecutive intersections with the section {a — 
0}. To avoid cumbersome notations related to the intersection of the domains of 
auxiliary maps, wherever appropriate we will assume that the Poincare maps are 
defined in a domain slightly larger than T> without further mention. Obviously 
there is no loss of generality in this assumption. 

We have seen in Proposition 17.31 that the Poincare map n preserves the 2-form 
A = G\(r, 9)rdr A d6, where 

G A (r,#) :=Bv a \ a=0 = l + O(e), 

and v a denotes the a-component of the local Beltrami field v. Mimicking the 
proof of this proposition, we immediately obtain that the Poincare map II of the 
divergence-free field u preserves the 2-form A := G.(r, 6)r dr A d9, with 

G A (r,fl) :=Bu a \ a=Q . 

Notice that u a does not vanish in {0} x V because the difference \u a — v a \ is small 
and the a-component of the local Beltrami field is close to 1. 
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Our next goal is to relate the above invariant 2-forms to apply Theorem 17.61 
More concretely, we will show that there is a C m -small diffeomorphism $ such that 

(7.30) A = $*A, 

where <&* is the pullback of the diffeomorphism. This will be done using Moser's 
trick. We start by noticing that the difference between these 2-forms is obviously 
exact, as 

A - A = dT 

with the 1-form T given by 

[G A {r,6)-G- A (f,6)]rdf \ d0 . 

Although we are making computations in polar coordinates, it is readily seen that 
all the objects we are considering are well-defined also at the origin, and therefore 
determine smooth forms in the whole disk V. 



Consider the non-autonomous vector field of class C°° 

jr{G A (f,8)~G A (f,9)}fdf 
s - r[(l-s)G A (r,0) + sG- A (r,e)} 



d r , 



where s will be the flow parameter. It should be noticed that, by the assumptions 
on the vector fields, the denominator behaves as r + 0(r 2 ) while the numerator is 
of order 0(r 2 ). The field Z s satisfies the C m bound 

(7.31) \\Z a \\c^(p) < CS 

for all s£ [0, 1] as a consequence of the estimate (I7.28[) . 

The time-s flow of the non-autonomous field Z s , which will be denoted by <j) a , is 
given by the solution to the initial value problem 

d 

—4> s x = Z s (4> s x), (j) x = x. 

OS 

Consider the s-dependent 2-form 

A s := (l-s)A+sA. 
A simple computation shows that 

(7-32) ^(#A.) =#(jU, + L 2 .A.) =<t>*(A-A + diz 3 As) = 0, 

where Lz s denotes the Lie derivative along Z s and the last equality follows imme- 
diately from the definition of A s and the fact that the interior product of Z s with 
A s is 

izAs = r. 

Therefore, if we set $ := 0i, we obtain Eq. (|7.30[) from (|7.32l) and the definition of 
A s . Moreover, 

(7.33) \\$-id\\ Cmiv) <CS' 

because <& is the time-1 flow of the vector field Z s , whose C m norm is controlled by 
Eq. (TOT]) . 

Let us now consider the map 

f[ := oflo$, 
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which is a diffcomorphism from a neighborhood of D 2 (which we still take as T>) 
onto its image. We shall next relate the new map II to the Poincare map II. By the 
definition of $ and the relation between the invariant 2-forms (|7.30[) . the map II 
preserves the same 2-form A as the Poincare map II, and II is close to II by (|7.33l) : 

(7.34) ||n-ft|| cm(1J) <CS'. 

We are now ready to apply Theorem 17.61 with the maps II and II. Indeed, the 
previous arguments and the way we have deformed the curve 7 ensure that the 
following statements hold true: 

(i) The rotation number cjn satisfies a Diophantine condition, so the map 
H\dj$2 is analytically conjugate to a rotation of frequency wn (cf. Proposi- 
tion EZ}. 

(ii) The normal torsion Afn is nonzero. 

(iii) The C k norm of II - ft is at most CS' by Eq. ([7T34]) . with 5' arbitrarily 
small. 

Hence Theorem 17.61 ensures that, for any given integer m, if the integer k is large 
enough and S' is sufficiently small there is a diffeomorphism ^ of T> with ^ — id 
arbitrarily small in the C m (T>) norm and supported in a neighborhood of 9D 2 such 
that the curve ^(c© 2 ) is invariant under the map II. Furthermore, the restriction 
of II to this invariant curve, n|^ 9D2 ), also has rotation number u)n. 

The definition of II then ensures that the curve $ o ^(<91D> 2 ) is invariant under 
the Poincare map II of the field u. (As an aside, notice that this invariant curve is 
close to <91D> 2 but, in principle, is not contained in the closure of D 2 , which is the 
reason why we are considering a slightly larger disk T> throughout the proof.) It is 
standard that this is equivalent to saying that there is an invariant torus of u whose 
intersection with the disk {a — 0} is precisely the aforementioned curve. Since the 
C m (T>) norm of the diffeomorphism $ o <f is arbitrarily small by the properties of ^ 
and Eq. (|7.33|) . one can take a diffeomorphism § of §} x V such that *(S] x 3D 2 ) is 
an invariant torus of the field u and — id is small in the C m norm and is supported 
in a neighborhood of §JxB 2 . Indeed, this diffeomorphism can be defined as follows. 
Take the solution to the system of ODEs 

dr u r (a,r,9) 
da u a (a, r, 9) ' 
d9 ug(a,r 7 9) 
da u a (a, r, 9) ' 

with the trajectory (r, 9) parametrized by the angle a and depending on the ini- 
tial conditions (r ,9g). Consider the function (p a mapping the initial conditions 
( r cb$o) to its time-a flow (r(a; ro, 9q), 9 (a; r , #o))- Then it is easy to check that 
the diffeomorphism ^ is given, in polar coordinates, by 

*(a, r, 9) := (a, ip a o $ o *(r, 6>)) . 

Actually, this formula simply asserts that the intersection of the invariant torus of 
the field u with each section {a — a Q }, understood as a curve in the disk T>, is the 
image under (p ao of the invariant curve $0 v[/(9ID) 2 ) at {a = 0}. Clearly the formula 
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for \I> and the estimates for the maps $ and '5 imply that 

II* ~ id||c™(s]xi>) < C\\$ ° * - id||c-(©) 
can be made arbitrarily small. 

To complete the proof of the theorem, it suffices to recall that the map n|^^ D2 j 
has rotation number u>u- This implies that the Poincare map II of the vector field u, 
restricted to the invariant curve $ o \P(9D 2 ), also has rotation number ton, which 
trivially implies that the vector field u itself is orbitally conjugate to a rotation of 
frequency oju on the invariant torus ^> (S) x <9B 2 ) . The existence of the diffeomor- 
phism of E 3 that appears in the statement of the theorem is then immediate. □ 

Remark 7.11. Ultimately, Theorem 17. 101 is a result on the preservation of invariant 
tori for small perturbations of the harmonic field h, rather than of the local Beltrami 
field v. Indeed, if we had considered a Beltrami field with small but otherwise 
arbitrary parameter A (possibly 0), we would have found the same expressions for 
the rotation number and normal torsion of its Poincare map, the only change being 
that the error terms would be C(e 2 + |A|) and 0(e 3 + |A|), respectively. 

We will conclude this section with a result on the persistence of another invariant 
set: we will show that, for an open and dense set of core curves 7 and small enough 
e, any divergence-free vector field that is a small perturbation of the local Beltrami 
field v has an elliptic periodic trajectory close to the core curve 7. We recall 
that a periodic trajectory is elliptic if the nontrivial eigenvalues of the associated 
monodromy matrix have all unit modulus but are different from 1. 

Proposition 7.12. Suppose that the total torsion of the curve 7 satisfies 

t(o) da nir 

for all integers n. Then for any S > 0, there exists some 8' > such that any 
divergence-free vector field u in the tube T e which is close to the local Beltrami 
field v in the sense that 

\\ u ~ v \\c k (S}xD^ < 6' • 

also has an elliptic periodic trajectory diffeomorphic to the curve {y = 0}. More- 
over, the corresponding diffeomorphism ^ is bounded by 

II* ~ id||c*($ixB 2 ) < s 
and is different from the identity only in a small neighborhood of the curve {y = 0}. 

Proof. From the expressions (|5.2[) - (|5.4[) for the harmonic field h, the estimates for 
the function tp proved in Theorem l5.1l and the connection between the local Beltrami 
field and h (see Eq. (|7.1|l ). we infer that the difference v — ho can be estimated as 

\\ v - ^o||c fc (SixB 2 ) < Cfe £ 1 

with the vector field ho defined as 

h := d a + t(o) (yi d 2 - 2/2 di) . 
It is clear that (a(s), y(s)) = (s, 0) is an ^-periodic trajectory of the field ho- Setting 

T := / r(a) da , 
Jo 
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an easy computation shows that the monodromy matrix of this trajectory for the 
field ho is 

/I \ 

cos T -sinT . 
\ sinT cos T J 

The nontrivial eigenvalues of this matrix are e ±lT , and hence different from ±1 by 
the hypotheses of the theorem, thus showing that {y = 0} is an elliptic trajectory 
of ho- 
lt is then standard that any divergence-free field u that is close enough to ho 
(say, || it — /io||c fc (sJxD 2 ) < <^i)) nas an elliptic periodic trajectory given by the image 
of {y = 0} under a diffeomorphism * with ||\I> — id|| C fc(§i xD 2) < S. There is no 
loss of generality in assuming that ^ — id is supported in a small neighborhood of 
{y = 0}. Since 

\\u - /lo||c fc (§ixB 2 ) < \\ u - v \\c k (S}xD 2 ) + \\ v - ^o||c fc (SixD 2 ) < + C/c e ! 

the theorem then follows by taking e and 5' small enough. □ 



8. Approximation by Beltrami fields with decay 

In this section we prove a result that allows us to approximate a field v that 
satisfies the Beltrami equation 

curl v = \v 

on a neighborhood of a compact set S, by a global Beltrami field u, which satisfies 

curl u = Xu 

in the whole space R 3 and falls off at infinity as 1/ 1 a? j . Throughout we will assume 
that the complement R 3 \S is a connected set. It is not hard to see that this 
condition is necessary. 

It will be more convenient for us to work with an auxiliary elliptic equation 
instead of considering the Beltrami equation directly. To this end, let us denote by 

cos\\x\ 
U[X) - Att\x\ 

the Green's function of the operator A + A 2 in ]R 3 , which satisfies the distributional 
equation 

AG + A 2 G = —So 

with So the Dirac measure supported at 0. We will use the notation Mr for the ball 
in R 3 centered at the origin and of radius R. 

The following lemma, which shows how to "sweep" the singularities of the 
Green's function, will be used in the demonstration of the global approximation 
theorem. Its proof is based on a duality argument and the Hahn-Banach theorem. 



Lemma 8.1. Take R > and consider a domain U C M 3 \!>2,r and a compact set 
S C Mr whose complement M?\S is connected. Let us consider the vector field 



M 



v ( x ) : = X! PmG(x - x m ) . 
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where {a; m }^_j is a finite set of points in Mr\S and p m £ R 3 are constant vectors. 
Then, for any S > 0, there is a finite set of points {zj}j =1 in the domain U and 
constant vectors Cj £ R 3 such that the finite linear combination 

J 

(8.1) w(x):=J2c j G(x-z j ) 

approximates the field v uniformly in S as 

\\v - w\\ c o(s) < S. 

Proof. Consider the space U of all vector fields that are linear combinations of the 
form (|8.ip . where the points Zj belong to the set U and the coefficients cj £ R 3 
are constant vectors. Restricting these fields to the set S, U can be regarded as a 
subspace of the Banach space C°(S, R 3 ) of continuous vector fields on S. 

By the Riesz-Markov theorem, the dual of C°(S, R 3 ) is the space 7W(5",R 3 ) of 
the finite vector-valued Borel measures on R 3 whose support is contained in the 
set S. Let us take any measure (X £ M(S, R 3 ) such that J R3 f-dp, = for all f £lA. 
Let us now define a field F £ Ll oc (R 3 , R 3 ) as 



F(x) := / G{x-x)dn(x), 
so that F satisfies the equation 

AF + X 2 F = -fi . 

Notice that F is identically zero on the open set U by the definition of the measure p, 
that R 3 \S f is connected and that F satisfies the elliptic equation 

AF + A 2 F = 

in R 3 \S'. Hence the unique continuation theorem ensures that the field F vanishes 
on the complement of S. It then follows that the measure p, also annihilates any 
field of the form p m G{x — x m ) because, as the points x m do not belong to S, 



= F(x m ) ■ p m = / G(X - X m ) p m ' dp(x) . 

Jr 3 

Therefore 

dp = , 



/ 



which implies that v can be uniformly approximated on S by elements of the sub- 
space U as a consequence of the Hahn-Banach theorem. The lemma then fol- 
lows. □ 

As an intermediate step before proving the global approximation result for the 
Beltrami equation, we will establish the following proposition on the approximation 
of solutions to the elliptic equation Aw = — \ 2 v by solutions defined in a large ball. 
Throughout, we will say that a differential equation holds in a closed set if it holds 
in a neighborhood of this set. 

Proposition 8.2. Let v be a vector field which satisfies the equation 
(8.2) Av = -\ 2 v 
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in a compact subset S o/R 3 . Assume that its complement M. 3 \S is connected and 
that S is contained in the ball Mr. Then for any 6 > and any positive integer k 
there is a vector field w satisfying the equation 

Aw = —\ 2 w 

in Mr that approximates the field v in S as 

(8.3) \\v - w\\ C k( S) < 6 . 
Here 6 is any fixed positive constant. 

Proof. By hypothesis, there is an open subset ft D S such that the field v satisfies 
the equation (I8.2j) in fi. We can assume that fi is contained in the ball Mr. Let us 
take a smooth function % : R 3 — > R equal to 1 in a closed set S' C fl whose interior 
contains S and identically zero outside tt. Defining a smooth extension v of the 
field v to K™ by setting v := x v i we obviously have 

(8.4) v{x) = [ G(x - x) p(x) dx 

with p := — Av — X 2 v. 

The vector field p is necessarily supported in tt\S' . Therefore, an easy continuity 
argument ensures that one can approximate the integral (I8.4[) in the compact set S' 
by a finite Riemann sum of the form 

M 

v(x) := ^ Pm G(x - x m ) 

m— 1 

so that, for any constant 5' > 0, 

\\v - i>\\c°(S') < 5' . 

Here p m are constant vectors in M 3 and x m are points that lie in Vl\S' . 

Let us take a domain U C R 3 \B2/j. Lemma [8.11 asserts that there is a vector 
field of the form 

J 

such that 

\\v - w\\ c °(s>) < 5' , 

where {zj}j =1 is a finite set of points in U and Cj £ R 3 are constant vectors. 
Therefore, 

(8.5) \\v-w\\ c o {sl) <26'. 

To complete the proof of the proposition, notice that the field v satisfies 

Av + X 2 v = 

in the set 5" (whose interior contains S) and w satisfies the same equation in the 
ball M2R. By standard elliptic estimates, it follows that the C° approximation (|8.5[) 
can be promoted to the C k bound 

h - w\\chs) < CkS' . 
Choosing 6' small enough, the result follows. □ 
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We are now ready to prove the global approximation theorem for the Beltrami 
equation with solutions that decay at infinity. To construct these solutions, we will 
truncate a suitable series representation for the fields in a large ball obtained using 
Proposition 18.21 and act on them using a convenient differential operator. 

Theorem 8.3. Let v be a vector field that satisfies the Beltrami equation 

curlu = Aw 

in a compact set S C K 3 , where X is a nonzero constant and the complement M. 3 \S 
is connected. Then there is a global Beltrami field u, satisfying the equation 

curl u = Xu 

in M 3 , which falls off at infinity as \D^u(x)\ < Cj/\x\ and approximates the field v 
in the C k norm as 

\\u-v\\ C h(s) < 8. 

Here 6 is any positive constant. 

Proof. Let us assume that the compact set S is contained in the ball Br/2- As the 
Beltrami field v satisfies the equation 

Av + X 2 v = 

in 5, by Proposition 18.21 there is a field w satisfying 

(8.6) Aw + A 2 w = 
in the ball Br and such that 

(8.7) ||w - w\\ C k+2 (s) < 5' . 

Let us take spherical coordinates (r, 9, tp) in the ball Br. Writing the field w as 
a series of spherical harmonics and using the equation (|8.6p we immediately obtain 
that w can be written as a series 

oo I 

w = ^2 ^ cimji(^r)Y lm (0,(p). 
1=0 m=-l 

Here ji is the spherical Bessel function, Yi m are the spherical harmonics and c; m € 
R 3 are constant vectors. Therefore, given any 8' > there is an integer L such that 
the finite sum 

l i 

^2 ci m ji (Xr) Yi m (0, ip) 

1=0 m=-l 

approximates the field w in L? sense: 

(8.8) \\u-w\\ L 2 {Mr) <6' . 

By the properties of spherical Bessel functions, the vector field u satisfies the equa- 
tion 

(8.9) Au + X 2 u = 
in M 3 and falls off at infinity as \D : 'u(x)\ < Cj/\x\. 

In view of Eqs. (|8.6I) and (|8.9|) . standard elliptic estimates allow us to pass from 
the L 2 bound ([8TS| to C k+2 estimate 

||u - H| cfc +2 (BiV2) < CkS' . 
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From this inequality and the bound (|8.7p we infer 
(8.10) ||fi- «||o*+»(s) < OS'. 

Let us now set 

u :- 



curl curl u + A curl u 



2A 2 

A simple computation shows that the vector field thus defined satisfies the Beltrami 
equation 

curl u = Xu 

in R 3 and falls off as < Cj/\x\ by the properties of u. Moreover, 

curl curl u + A curl u 



\\u — v\ 



C k (S) 



2A 2 

(curl+A) curl(u — v) 



C k (S) 



<C\\u- 
< C5' 



2A 2 

v \\C k + 2 {S) 



C k (S) 



as we wanted to prove. 



□ 



Remark 8.4. The fall-off at infinity of the global Beltrami field u is obtained from the 
truncation of the explicit series representation for the auxiliary field w. This is the 
reason why Theorem 18 .31 does not work in arbitrary open Riemannian 3-manifolds, 
unlike the approximation theorem we used in |llj . Notice that the latter theorem 
does not yield any control at infinity whatsoever for the global Beltrami fields. 



9. Proof of the main theorem 



We are now ready to give the proof of Theorem ll.il Let us begin by considering 
one of the curves (say, 71) in the statement of the theorem. By perturbing this 
curve with a C m -small diffeomorphism if necessary, we can assume that 71 is an 
analytic curve whose curvature k does not vanish anywhere p. 184], and that its 
total torsion satisfies 

i 

t(ol) da ^ 71" k 

for all integers k. As before, I denotes the length of the curve 71. It is worth 
emphasizing that, as these conditions are open, they are obviously preserved if we 
deform the curve 71 with a diffeomorphism that is close enough to the identity in 
the C m norm (m ^ 3). 

Let us consider the tube % (71) of core curve 71 and thickness e, and the corre- 
sponding harmonic field h, given by the expression (I5.3[) in the coordinates (a, y) 
adapted to the tube. Throughout we will assume that e is small enough. By The- 
orem [6751 we can consider the associated local Beltrami field i>i, which is given by 
the only solution to the Beltrami equation with parameter A = e 3 , 

curlvi = e 3 vi , 

in the tube 7e(7i) that is tangent to the boundary and whose harmonic part is 
the field h. (Notice that in this section Vi will not stand for the component of a 
vector field in the direction of the coordinate y\.) Since the boundary of the tube 
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is analytic, it is well known l 23: that the local Beltrami field V\ is analytic in the 
closure of a small neighborhood Q,\ of 71 (71). 

Since the total torsion of the curve is not an integral multiple of 7r, Theorem 1 7. 101 
and Proposition 17. 12l ensure that, given any S > 0, we can deform the curve 71 by a 
diffeomorphism of M 3 arbitrarily close to the identity in the C m norm so that any 
divergence-free vector field u in VL\ with 

II - "- u l|lc fe (Oi) < s ' 

has: 

• An invariant tube given by ^1 [7^(71)]. 

• An elliptic periodic trajectory given by ^1(71). 

Moreover, on the invariant torus d^>i[%(ji)], the field u is orbitally conjugate to a 
Diophantine rotation, and therefore ergodic. Here ^1 and V^i are diffeomorphisms 
of R 3 with 

11*1 - id||c m (R 3 ) + ||*i - id||c m (K 3 ) < 5 
and such that the differences "fi — id and *i — id are supported on small neighbor- 
hoods U\, U\ of 97^(7i) and 71, respectively. The constants k and 5' depend on m, 
S and on the geometry of the curve 71 . 

We can apply the same argument for each curve 7$ (1 ^ i ^ N), thereby ob- 
taining (for small enough e) a collection of local Beltrami fields Vi, satisfying the 
equation 

curlui = e 3 Vi 

in the closure of a neighborhood f2, of the closed tube TJcii) and such that any 
divergence-free vector field u in Qi with \\u— t'i || c fe (£ii) < has an invariant torus, 
where the field u is ergodic, and an elliptic periodic trajectory. Furthermore, they 
are respectively given by d^i[T e (ji)] and ^(7.;), where ^>i, are C m -small diffeo- 
morphisms of M 3 with 'J'i — id and — id supported in small neighborhoods Ui, Ui 
of dTe^i) and 7$, in each case. We can assume that the complement 

M 3 \(f7iU---Uft7v) 
is connected and that the sets fit are pairwise disjoint. 

Let us define a vector field v in f2i U • • • U fijv by setting it equal to the local 
Beltrami field Vi in each set 0,-. By Theorem 18.31 there is a Beltrami field u, which 
satisfies the equation 

curlu = ( 3 u 

in R 3 , that falls off at infinity as \D^u(x)\ < Cj/\x\ and approximates the field v as 

II" ~ «llc fc (^iu-uf2„) < ■ 
Therefore, if we define the diffeomorphism $ of R 3 as 

(%(x) ifxeUi, 
%(x) ifxe&i, 
x otherwise , 

it follows that, for each i, $[7I(7i)] is a vortex tube of the Beltrami field it, and that 
<I>(7i) is an elliptic periodic trajectory. Besides, the Beltrami field u is ergodic (and 
orbitally conjugate to a Diophantine rotation) on each invariant torus <9<I>[7^(7i)]. 
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The field u being orbitally conjugate to a Diophantine rotation on each invariant 
torus <9$[7^(7i)]) it follows [14] that this invariant torus is accumulated by a Cantor- 
like set of invariant tori with positive Lebesgue measure. On these invariant tori, 
the field is also orbitally conjugate to Diophantine rotations. The corresponding set 
of Diophantine frequencies is Cantor-like because the normal torsion An is nonzero. 
Therefore, if we consider the trajectories of the field u between two of these invariant 
tori (lying on the same vortex tube $[7^ (7t)]), Angenent's dichotomy [1, asserts that 
either there is a horseshoe-type invariant set between them or there are invariant 
tori where the field is conjugate to a rotation of rational frequency. In both cases, 
there are infinitely many periodic trajectories between these tori. The theorem 
then follows. 

Remark 9.1. It is worth giving some additional details about what we understand 
by a thin tube. What we have proved is that for any set of smooth periodic curves 
7i, . . . ,7jv, S > and meN, there is a positive constant eo, depending on S, m 
and the geometry of the curves, such that the statement of Theorem 1 1 . 1 1 holds true 
for any thickness e < eo, the diffeomorphism $ that maps the tubes 7l("fi) into 
vortex tubes being bounded by 11$ — id||(7m( E 3) < i5. It should be noticed that both 
Theorem ll.il and its proof hold verbatim if one takes tubes 7^(71), . . . ,% N ("fN) of 
different thickness, as long as ei, . . . , ejv are small enough (that is, smaller than the 
above constant eo). 

Remark 9.2. Contrary to what happened in the main theorem of [11) . the proof 
of Theorem 11.11 does not work in general (even if we drop the requirement that 
the field u decays at infinity, replacing Theorem 18.31 by [TT| Theorem 3.6]) if we 
substitute the finite set of curves {%}fLi by an infinite set {7^}^ that is locally 
finite. The reason is that the "maximal thickness" £0(7*) associated to each curve 7* 
individually does not need to be bounded away from zero: since all the vector 
fields Vi, defined in a neighborhood of the tube T ei {ji) (whose thickness can vary 
from tube to tube), must satisfy the Beltrami equation 

curium = Xvi 

with the same constant A for all i. Since this A must be of order 0[eo(7i) 3 ] for all i, 
it is clear that the construction breaks down if the infimum of the positive quantities 
{ e o{ji)}iZi is 0- On the other hand, if this infimum is positive, we can apply the 
approximation theorem in [11] Theorem 3.6] to construct a global Beltrami field 
with all the properties listed in the statement of Theorem 11.11 with the exception 
that its growth at infinity is not controlled. 

10. A COMMENT ABOUT THE NAVIER-STOKES EQUATION 

To conclude, we will present an easy application of Theorem II .11 to the existence 
of (time-dependent) solutions to the Navier-Stokes equation that have a prescribed 
set of stationary (possibly knotted and linked) vortex tubes. 

For this, let us take the global Beltrami field u that we considered in Section |9l 
That is, u satisfies the equation 

curlu = e u 

in R 3 , falls off at infinity as C/\x\ and has a set of thin invariant tubes given by 
{^[Te( r yi)]}iLi, where the diffeomorphism 3> is close to the identity and e is a small 
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constant. Then the analytic time-dependent field 

w(x,t) := u(x)e^^ 6 * 

is a solution of the Navier-Stokes equation 

dtW + (w • V)ui = vAw — VP , div w = 

in M 3 with pressure P(x,t) — c — ^\w(x, t)\ 2 . As the vortex lines of w (which are 
the trajectories of the vorticity curl w(x, t) for fixed t) coincide with those of u at all 
times, up to a reparametrization, it follows that w is a solution to the Navier-Stokes 
equation with the desired properties. 
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